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ABSTRACT 

An isoparametric finite element formulation solving a 
general form of the transient field equation is presented in 
this thesis. The formulation is developed for fields in 
three-dimensional Euclidean space. A FORTRAN IV computer 
program employing double precision arithmetic, compact 
storage techniques, and providing the option of several 
numerical integration methods and types of finite elements 
is presented. Data input may be in rectangular, cylindrical, 
or spherical coordinates and variations in time integration 
step size are permitted. Time dependent boundary conditions 
are not considered. 

Comparisons of theoretical and computer solutions for a 
variety of test problems demonstrate close agreement. 
Instructions for use of the program are included. 
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I. 



INTRODUCTION 



The 'quasi-harmonic' partial differential equation: 






+ 



9z 



k 



z 9z 



3d> 9 2 t}> 



= 0 



( 1 ) 



applies to a variety of physical problems [Ref. 10]. The 
Laplace equation is a typical and well known example of a 
specialized form of this equation. Since an analytical 
solution in closed form is seldom available for useful 
engineering problems, the development of a practical general 
purpose computerized solution method represents a useful 
addition to the analytical tools of the engineer. 

This thesis describes a general purpose FORTRAN IV 
computer program using the isoparametric finite element 
approach to solve the field equation in general form. The 
program does not presently permit a time variation of boun- 
dary conditions or a functional variation of material 
properties . 

The program allows selection of the linear, quadratic, 
or cubic isoparametric finite element with the option of 
specifying two through six Gauss points for numerical inte- 
gration in space. Data input may be a rectangular, cylin- 
drical, or spherical coordinates. Five time integration 
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routines are provided. Output is presented in a tabular 
form of node number, space coordinate, and functional value. 

The structure of the program is modular in order to 
accommodate improvements and adjustments due to new 
developments . 

This program was constructed and tested on an IBM-360/67 
computer system with OS/360 release 18. Other systems and 
installations will require modifications. 
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II. FINITE ELEMENT FORMULATION 



The 'quasi-harmonic' or field equation describing the 
behaviour of some physical quantity, 6, is: 



_3_ 

3x 





+ 




(1) 



Commonly imposed boundary conditions are of the forms: 



4> = 4> b 



( 2 ) 



k 2 + 

x 3x x 



v y 9 y 



ii a 



+ k 



f± i 

3 z : 



+ q + a<f> = 0 



(3) 



Where <J> D represents the value of 4> specified on the boundary; 

Jj 

J 1,1, and Z are the direction cosines of the outward 
x y z 

normal to the boundary surface; q is an outward flux per 

unit area; Q a generation per unit volume; and k , k , k , 

x y z 

y, and p are material properties. Q, k , k , k , y, and p 

x y z 

may be functions of space and time. 

Reference 10 thoroughly develops the concept of a finite 
element and the process leading to a discretized form of 
equation (1) and the boundary conditions given by equations 
(2) and (3). If the scalar variable, <J) , is described by a 
shape function, <N>, and nodal values of the variable, {<}>}, 
such that: 



<j> = <N> {<J>} 



(4) 



the result is: 



[H] { } + [C] { 4>} + [ G ] { 4> } + {F} = 0 (5) 

[H] , [C] , and [G] are symmetric matrices defined, on the 

element level, as: 

[C] e = / <N> T y <N> dx dy dz (6) 

e 



[G] e = / <N> T p <N> dx dy dz 



(7) 



r „,e r f u „ 9 *k . u 

[H] = / k <v— > <~— > + k <~— > > 

1 v \ x 3 x 3 y 3 3y 

e \ x J y 2 



. , „3N* 3N. \ , , , 

+ k z <— : > < ~ > J 



( 8 ) 



The load vector describing boundary conditions is: 

{F} e = -/ Q<N> T dxdydz + / q <N> T ds 
v e e 

+ <N> T a <N> ds ^ {<f>} e (9) 

where the integrals are performed only over the surface or 
volume where the prescribed boundary condition applies. 

Element matrices and the load vector are assembled by 
standard methods to form equation (5) which is the discre- 
tized form of equation (1) . Adding the load vectors given 
by equation (9) to the rest of the nodal intensities 
guarantee the satisfaction of the boundary conditions given 
by equation (3) . Finally, the boundary condition of equa- 
tion (2) can be applied to the assembled set of equations. 
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The finite element technique models a continuum as an 
assembly of substructures; thus, for an element of small 
size, the parameters k^, k , k^, jj , p, a, Q, and q may be 
considered constant and the integrations required by equa- 
tions (6), (7), (8), and (9) reduced to an automated numeri- 

cal technique. 

A. ISOPARAMETRIC FINITE ELEMENTS 

The linear, quadratic, and cubic isoparametric shape 
functions are utilized in this development. The resulting 
elements contain 8, 20, and 32 nodal points respectively. 

A full development of the concept and properties of this 
family of finite elements is contained in Ref. 10; only the 
necessary results will be repeated here. 

A "local" non-dimensional system of reference (C,n,?) is 
used for convenience of calculation. Results are transformed 
into the "global" cartesian (x,y,z) system of reference prior 
to final assembly of problem matrices and vectors. The 
shape functions, N^(£,n,?)/ for the element types used in 
this development are given in Ref. 10 and need not be 
repeated here. The nodal numbering conventions employed for 
each element type are shown in Figures 1, 2, and 3. 

Important results relating the shape function to the 
local and global reference systems [Ref. 10] are: 
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Figure 1. Linear Isoparametric Flement. 

5 




Figure 2. Quadratic Isoparametric Element. 
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Fiaure 3. 



Cubic Isoparametric Element. 
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where the Jacobian matrix, [J], is defined as 
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and 



( 10 ) 



( 11 ) 



dxdydz = det [J] d£ dri d£ (12) 

The partial derivatives of the shape function with 
respect to global coordinates are recovered from equation 
(10) by applying the inverse of the Jacobian matrix. 



B. MATRIX FORMATION 

Assuming constant material properties within an element 
and applying equations (10), (11), and (12) to the matrix 

and vector formulations of equations (6), (7), (8) and (9) 

reduces the form of the integrals to: 

/ <n> T <N> det [J] d£dnd£ 

$ 

T 

f <|^> <|^> det [J] dCdr|d;+... 

f <N> T det [J] d£dr|dc 
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/ <N> T det [J ] ds 
s e 

/ <N>^ <N> det [J ] ds 

s e 

where the integrations are performed in the local nondimen- 
sional coordinate system. The forms of equations (6) , (7) , 

(8), and (9) become ideally suited for a Gaussian numerical 
integration procedure. 

C. BOUNDARY CONDITIONS 

The specification of a value of the function ({> along a 
given boundary as in equation (2) is straightforward and 
requires no comment. 

Boundary specifications of equation (3) appear in the 
expression for a forcing vector given by equation (9). The 
first term of equation C9 ) containing Q represents a gener- 
ation rate per unit volume; the second term containing q 
represents a flux per unit area; and the third term contain- 
ing a represents a surface loss coefficient per unit area. 
The terms containing q and Q contribute to a forcing vector, 
{F}, while the term containing a results in a correction to 
the corresponding element [H] matrix. 

If, in equation (5) , q and a are zero and k , k , and k 

x y z 

are equal the boundary condition reduces to the requirement 
that the normal derivative be equal to zero. This condition 
is automatically imposed by equation (5) in the absence of 
any other boundary condition specif ication . 

Arpaci [Ref. 1] gives a number of boundary conditions 
applicable to the heat problem. Of these the Radiation, 
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Interface in Motion, and Change of Phase are not included in 
this development. 

Consider an element located at the surface of a body in 
space. From Figures 1, 2, and 3 observe that a corner node 
may be common to one, two, or three boundary surfaces per 
element and that a mid side node is common to one or two 
boundary surfaces per element. Each boundary surface may 
be subject to a different boundary condition specification. 

If two or more of the element surfaces have different boundary 
conditions imposed, a unique situation arises if the boundary 
specifications include conditions on 4 > and q or a. In this 
case the boundary conditions must be specified in order. 
First, all boundary conditions involving q and ct are formed 
and merged into the overall system of equations; then the 
boundary conditions involving 4> are applied. 

D. SOLUTION TECHNIQUES 

For a fixed point in time and for steady state problems, 
equation (5) reduces to the form 

[A] { 4> } = {B} (13) 

where [A] is a n x n square matrix. The vector (b) is fully 
determined in all cases except where a nodal value of 4> is 
specified. If a nodal value, 4 k, is specified the corres- 
ponding forcing vector term, , will contain the unknown 
quantity of the equation. 
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1. Systems of Linear Equations 

Equation (13) may be solved by standard methods of 
linear algebra. The method employed here to sort the known 
and unknown quantities of equation (13) was developed by 
Irons [Ref. 6] and has the advantages of preserving the 
symmetry of matrix [A] and permitting the determination of 
all unknown quantities of the system of equations. 

After applying the Irons [Ref. 6] modifications, 
solution of equation (13) is accomplished by the reduction 
of matrix [A] into the product of lower unit triangular, 
diagonal, and upper triangular matrices such that: 

[A] = [L] [D] [L) T 

Unknown quantities are found by performing a forward sub- 
stitution followed by a backward substitution [Ref. 13). 

2 . Steady State 

In steady state, equation (5) readily reduces to 
the form of equation (13) and is solved by the technique 
described previously. 

3 . First Order in Time 

At fixed points in time, equation (5) reduces to the 
form of equation (13) . The numerical integration schemes 
employed are straightforward and are summarized below. The 
initial value of {<}>} is assumed known and time step size is 
denoted by h. 

a. Euler I Method 

This integration scheme is well documented by 
Crandall [Ref. 3). The result, for the i*"* 1 point in time, is 
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[C] { 4»} 1_1 = -[Hit.})} 1 " 1 - { F} 


(15a) 


U) 1 = { 4>) 1-1 + h U) 1-1 


(15b) 



by Ref. 



b. Trapezoidal Rule Method 

This approach to time integration is documented 
3 and Ref. 11. At mid interval 

{4,} i+3 s = U) 1 + U) i+1 



h 



When applied to equation (5) , the resulting recursion 
formula is: 

( [H] + i [C]){<})} 1+1 = (- [H] + ^ [C]){(j)} 1 - 2{F} (16) 

h h 

c. Linear Time Shape Function Method 

Reference 10 develops the following approach to 
time integration. Assume that {<{>} may be considered as a 
product of the nodal values of {<t>} at discrete times and a 
time shape function, M(t) , such that 

<M> = <M q ( t) ,M (t) > 
where, for a linear variation: 



and 

M x = t/h 
the result is: 
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{ 4> } = <M , M x > 



uj 1 

Applying this relationship and the Galerkin process described 
in Ref. 10 results in the formula: 

(| [H] + i [C] ) U } 1 = - (3 [H3 - i [C] ) { 4>) 1-1 - {F} (17) 

4 . Second Order in Time 

At fixed points in time, equation (5) reduces to the 
form of equation (13) . The integration methods employed are 
summarized below. Initial values of {<}>} and {<J>} are assumed 
to be known . Time step size is denoted by h. 

a. Euler I Method 

The Euler I integration method is well known and 
documented by many standard sources such as Ref. 3. The 
result, when applied to equation (5) is: 

[G] (i ) 1-1 = - [HI U ) 1 " 1 - [CIU ) 1-1 - {F} (18a) 

{ 4> > 1 = {<}>} 1-1 + h {4>} 1_1 . (18b) 

U} 1 = -U} 1 " 1 + h {cj>} 1-1 (18c) 

In the event that the [C] matrix equals zero these equations 

remain valid. 

b. Finite Difference Method 

Finite difference approximations are developed 
in Ref. 3. A four point backward difference approximation 
for {<}>} and a three point backward difference approximation 
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for {<}>} have truncation errors of the same order of magni- 
2 

tude, h , and are: 

(i) 1 = (-{4>} 1-3 + 4{ 4>> X_2 - 5{4>} 1-1 + 2{<t>} 1 ) (19a) 

h 

{ ^ }i = K {( ^ }i " 2 - 2{^} i_1 + § U) 1 ) (19b) 

Define 

(4>* } = (-{ <})} i_3 + 4 { <J> } i— 2 - 5{4>} i " 1 ) (19c) 

h 

i 

and 

{<*>**} = £ (| U) 1 " 2 - 2{ <}>} i_1 ) ( 19d) 

Equations (19a) and (19b) when applied to equation (5) give 
the following recursion formula 

( [H] + \ [Cl + [C] ) { 4>> 1 = 

h 

- {F} - [G] {<()*} - [C] {<*>**} (20) 

if [C] is not equal to zero and 

( [H] + [G]) { 4>> 1 = - {F} - [G]{<|>*} (21) 

h 

if [C] equals zero. 

Equation (21) is known as the Houbolt method and is 
shown by Johnson [Ref. 7] to be unconditionally stable. 

Three starting values of {<{>} must be obtained by inde- 
pendent means such as the application of the Euler I method 
or a Taylor series expansion. 
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III. DESCRIPTION OF THE PROGRAM 



A functional block diagram of the program is presented 
in Figures 4 through 9 . A modular structure was adopted in 
order to permit simple modifications to accommodate future 
improvements and the adaptation of the program for special 
purposes . 

The following paragraphs give a brief description of the 
function and operation of the program. Appendices A through 
D discuss the preparation of input data, list the important 
computer program nomenclature, and list the program. A full 
understanding of the operation of this program can only be 
achieved by a study of Appendices A through D and an under- 
standing of the principles and solution techniques employed. 

A. ASSUMPTIONS AND RESTRICTIONS 

The key assumptions and restrictions which impose 
restrictions in the type of problems that can be solved by 
the program are: 

1. Material properties k^, k^, k z , p, U, and a and 
problem parameters q and Q are constant functions of time 
and space throughout the element and are independent of the 
function <f>. 

2. Boundary values of $ are not functions of time. 

3. The principal axes of each element align with the 
principal axes of the global coordinate system and are fixed 
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MAIN 




Figure 4. Subroutine Calling Sequence for 
Steady State Problem. 




Figure 5. Subroutine Calling Sequence for 
First Order Time Problem. 
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NOTE: MERGE (2) is not called if the [C] 

matrix is absent 



Figure 6. Subroutine Calling Sequence for 
Second Order Time Problem. 
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NOTE: STUFF not called for the Steady State Problem 



Figure 7. Assembly Routines for the Formation 
of Matrix [H] . 
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STUFF only called 
from MERGE (2) for 
the Second Order 
Time Problem 



Figure 8. Assembly Routines for the Formation 
of Matrices [C] and [G] . 




NOTE: (1) STUFF called by TIME23 only 

(2) MULT not called by STEADY 



Figure 9 . Service Routines Used by 
Solution Routines. 
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in direction throughout the element. 



These limitations will be discussed in Section V. 

B. GENERAL PROGRAM ORGANIZATION 

The program may be divided into the three major indepen- 
dent functions: reading of input data, matrix formation, 

and actual problem solution. This organization is shown in 
Figures 4 through 9. 

All matrices involved in the problem are symmetrical. 

At the element level the complete matrix is formed; however, 
formation of the master matrices [II] , [C] , and [G] and prob- 

lem solution utilize a standard compressed storage mode to 
conserve computer core. The technique consists of placing 
the diagonal elements of each row of the matrix into the 
first column and retaining only the upper elements of the 
matrix . 

Reading of necessary input data is accomplished by sub- 
routine INPUT. Preparation of input data is discussed in 
Appendix B. An echo check of all input data is provided. 

Formation of the required matrices and vectors is 
controlled by subroutine MERGE. Matrices [H] , [C] , and [G] 

are formed sequentially. This permits output of the 
completed matrix to an external storage device and conserves 
core storage. Boundary condition calculations are performed 
during formation of the [H] matrix. Associated subroutines 
include CUBE, FORM, and JACOB. 
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The transient solution subroutines are similar in con- 



struction and differ only in the integration procedure 
employed. Time integration step size changes and data 
output are controlled by the solution subroutines. Certain 
functions common to all solution subroutines are performed 
by a set of service subroutines which include LDLT, BNSN, 
STUFF, PUNC, and MULT. 

C. SUBROUTINES 

The MAIN program is simply a control segment that serves 
to call INPUT; MERGE(l), MERGE ( 2 ) , or ME RGE ( 3 ) in the proper 
order and as required; and the applicable solution subroutine 
STEADY, TIME 11 , TIME13, TIME21, or TIME23. No calculations 
are performed in the MAIN program 

1 . Subroutine INPUT 

This subroutine reads all input data and provides the 
appropriate echo check. The number of nodes per element is 
determined from the element type specification and, if 
required, cylindrical or spherical coordinates are trans- 
formed into rectangular coordinates by calling subroutine 
MARK. 

2. Subroutine MERGE (IMTRX) 

This subroutine forms the master matrix specified 
by the calling argument, IMTRX, by assembling each element 
matrix into the master matrix. A call to subroutine CUBE 
returns the desired element matrix. IMTRX=1 corresponds to 
the [H] matrix, IMTRX=2 corresponds to the [C] matrix, and 
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IHTRX^S corresponds to the [G] matrix. Subroutine BCOND is 
called during formation of the [H] matrix and, if required 
by the solution subroutine selected, the [H] or [C] matrices 
are placed into the appropriate storage location after 
formation . 

3 . Subroutine BCOND 

This subroutine forms the element vectors and matrices 
required by the appropriate boundary condition specified and 
assembles the resulting vector or modification matrix into 
the forcing vector, { F } , or the corresponding [H] matrix. 
Formation at the element level is accomplished by a call to 
the appropriate entry point of subroutine CUBE. 

4 . Subroutine STEADY 

This subroutine solves the steady state problem. 

After the application of appropriate nodal boundary condi- 
tions, a call to subroutine LDLT and entry point SLV produces 
the solution. Data output is a printout of the (4>) and { F} 
vectors and punched card output of the {$ vector if requested. 

5 . Subroutine TIME11 

This subroutine employs the Euler I time integration 
method, equations (15a) and (15b), to solve the first order 
time problem. One call to subroutine LDLT is required; sub- 
sequent calls to entry point SLV produce the desired solu- 
tion. Data output is a tabular presentation of the 
vector at specified points in time with punched card output 
of the {<t>) vector provided if requested. 
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6. Subroutine TIME13 



This subroutine employs the "Linear Time Shape 
Function" as the time integration method (KINT=3) , equation 
(17) , or the Trapezoidal Rule time integration method 
(KINT=2) , equation (16) , to solve the first order time prob- 
lem. External core storage devices are utilized to store 
the [H] and [C] matrices in order to permit restarting the 
problem at each change in time integration step size. One 
call to subroutine LDLT is required for each change in time 
integration step size with the required solution obtained by 
sequential calls to entry point SLV. Data output is the 
same as for TIME11. 

7 . Subroutine TIME21 

This subroutine employs the Euler I time integration 
method, equations (18a), (18b), and (18c), to solve the 

second order time problem. Operation of this subroutine is 
identical to the operation of TIME 11 except for the inclusion 
of the vector {$} in the data output. 

8 . Subroutine TIME23 

This subroutine employs a Four Point Backward Finite 
Difference (Houbolt) time integration method, equations (20) 
and (21) , to solve the second order time problem. Starting 
values are obtained by a five term Taylor Series expansion. 
Operation of this subroutine is similar to the operation of 
subroutine TIME13 employing external core storage devices 
and saving the values of {<f>} required to restart the problem 
after a change in time integration step size. Data output 
is the same as for TIME21. 
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9. Subroutine CUBE (NPEL , NGP , ITYPE ) 



This subroutine performs the Gaussian numerical inte- 
gration required for the formation of the various element 
matrices. Entry arguments are the number of nodes per 
element, number of Gauss points specified for the integration 
process, and the type of integration to be performed. Gauss 
points and weighting factors were obtained from Ref. 9. A 
selection of two through six Gauss points with a default 
option of three Gauss points is provided. For each specified 
Gaussian integration point, the numerical evaluation of the 
integrand is determined by calling subroutine FORMC at the 
appropriate entry point and applying the proper weight 
factor to the result prior to summation. 

The specialized volume and surface integrals required 
for boundary condition calculations are obtained by entering 
this subroutine at entry point CUBEB. 

10. Subroutine FORMC (XI , ETA, ZETA,NPEL) 

T 

This subroutine performs the products of <N)>, 

,3M x T.8Nv . , , . . . , , 

y ' etc. required at each point in local coor- 
dinate space for the numerical integration process. Values 
of the shape function and partial derivatives of the shape 
function with respect to the local coordinate system are 
determined by a call to subroutine SHAPE. Calling subroutine 
JACOB returns the required values of [J], [J] ^ , and Det[J]. 

Necessary entry arguments are the local coordinates of the 
point in space under consideration and the number of points 
per element. Various entry points are employed to accommodate 
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the specialized calculations required in the formation of 
each element matrix or vector and volume or area integrations. 



11 • Subroutine SHAPE (XI .ETA, ZETA ,NPEL) 

This subroutine forms the shape function and the 
partial derivative of the shape function with respect to the 
local coordinate system at the point in local coordinate 
space specified by the entry argument. The type of element 
to be considered is determined by the specification of the 
number of points per element in the calling argument. 

12. Subroutine JACOB (NPEL , AJ , DCOL , CORD, AJI , DTJ) 

This subroutine evaluates the Jacobian matrix, 

determines the inverse of the Jacobian matrix and the deter- 
minant of the Jacobian. The number of points per element, 
global coordinates of the node points, and shape function 
partial derivatives with respect to the local system of 
reference are required entry arguments. [J] is returned in 
AJ, [J] in AJI, and det[J] in DTJ. 

The specialized calculations required for a surface 
integration are accomplished by entering this subroutine at 
entry point JACOBB. 

13. Subroutine LDLT (A ,N ,M, ASN) 

This subroutine reduces the matrix [A] into the 

products of lower unit triangular, diagonal, and upper tri- 

T 

angular matrices, [L] [D] [L] , required by the process of 

solving the linear system of equations described by 
[A] {x} = {B}. N and M are the number of equations and half 

band width respectively. ASN is an arbitrarily small number; 
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numbers less than ASN squared are considered to be zero. 

Entry point SLV(Z,B,N,M) performs the forward substitution 
and backward substitution required to solve the system of 
equations . 

14 . Subroutine MARK(C,K,M) 

This subroutine converts coordinates given in a 
cylindrical or spherical system of reference to a rectangular 
system of reference. C is a matrix containing the specified 
coordinates, K is a code to indicate the system of reference 
employed, and M is the number of coordinates specified. 

15. Subroutine BNSN (A , N , APZRO , ABIGN) 

This subroutine determines the arbitrarily small 
number required by subroutine LDLT and the arbitrarily large 
number required by the equation solving process. APZRO is 
determined from the matrix [A] by finding the average value 
of the diagonal elements and dividing the result by 10^ 3 . 

ABIGN is 10 30 x ASN. This provides the numerical spread 
necessary for the IBM/360 computer installation employed for 
this program development. 

16. Subroutine MULT (A , B , C , N , M) 

This subroutine performs the matrix multiplication 
[A] {b} = {C} in a compressed storage form. This multipli- 

cation is required by the various time integration subroutines. 
N and M are the number of equations and half band width 
respectively . 
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17. Subroutine STUFF (A , B , N , M) 

This subroutine places the matrix [A] into matrix 
[B] ; or entry STUFFV (AA , BB , N , MC) places the vector (b) into 
the specified column, MC , of matrix [AA] . 

18 • Subroutine PUNC (A , B , N , K) 

This subroutine produces the punched output requested 
by the solution subroutines STEADY, TIME11, TIME13, TII1E21, 
or TIME 2 3 . 



I 
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IV. 



TEST PROBLEMS 



The test problems employed here are designed to test and 
verify the operation of the computer program. The four cate- 
gories of tests performed were assembly, steady state, first 
order in time, and second order in time problems. 

The assembly tests served to verify the calculations and 
integrations required to form the problem matrices and vectors 
on the element level. Steady state problems were used to 
verify the matrix assembly process, linear equation solution 
process, and the validity of boundary condition formulations. 
Transient problems of first and second order in time served 
to test the various solution subroutines utilized. 

A. ASSEMBLY TESTS 

A rectangular region in space provides a convenient means 
of verifying element matrix and vector assembly procedures . 

The calculation of all problem matrices and vectors is 
straightforward and element calculations of the [H] , [C] , and 

[G] matrices were verified in this manner. Similarly, the 
basic assembly processes of subroutines MERGE and BCOND are 
readily verified. 

The integrations required by equations (6) , (7) , (8) and 

(9) are performed by a Gaussian numerical integration scheme. 
In the Gauss procedure, n sampling points will exactly 
integrate a polynomial of degree (2n-l) . The integrations 
necessary may involve polynomials of degree 2 through 14 in 
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my given direction of the local coordinate system. Poly- 
lomial degree will depend on the shape of the clement 
elected. For practical applications, a selection of two 
hrough six Gauss points is sufficient. For a rectangular 
hape the .linear quadratic, and cubic elements are exactly 
ntegrated by 2 , 3, and 4 Gauss points respectively. Veri- 

fication of the number of Gauss points required to exactly 
.ntegratc elements of other shapes was not considered. 

5. STEADY STATE PROBLEMS 

1 . Axial Temperature Distribution in a Cylinder 

Consider a cylinder of the geometry specified in 
^igure 10. If the sides arc perfectly insulated a linear 
semperature distribution results when end conditions arc 
specified [Refs. 1 and 5]. Specify the thermal conductivity, 
as 20 BTU/hr/ft/°F and divide the region into two finite 
ilements. The top surface is subjected to an ambient tem- 
perature of 1 0 0 ° F . 

If the bottom of the cylinder is subjected to a 

2 

ieat flux, q, of 500 BTU/hr/ft or maintained at 300°F with 

2 

i heat transfer coefficient, a, of 5 BTU/hr/ft /°F identical 
icmperaturc distributions result. Theoretical and finite 
ilement solutions are shown in Table I. 
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Figure 10. Axial Cylinder Problem Geometry. 



TABLE I 

AXIAL TEMPERATURE DISTRIBUTION IN A CYLINDER 



Distance 

from 

Bottom (ft) 




Temperature (°F) 




Theoretical 

Solution 


2 Quadratic 
Element 
Mode 1 


2 Cubic 
Element 
Model 


0.00 


200 . 00 


200.00 


200 . 00 


0.67 


183. 33 




183.33 


1.00 


175.00 


175.00 




1.33 


166.67 




166.67 


2.00 


150.00 


150 . 00 


150.00 


2.67 


133.33 




133.33 


3.00 


125.00 


125.00 




3.33 


116.67 




116.67 


4.00 


100.00 


100.00 


100.00 
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2 . Radial Temperature Distribution in a Cylinder 

Consider a hollow cylinder of infinite length. A 
segment of such a cylinder is shown in Figure 11. Heat is 



(lost by convection through the outer wall. Consider the case 
of the interior surface maintained at a constant temperature 
or the case of a specified heat generation within the cylinder 
rfith the interior surface perfectly insulated. Assumed 
problem parameters are: 



a . 


Inside radius (r 


3 in . 


b. 


Outside radius (r Q ) 


9 in . 


c. 


Heat transfer coefficient (a) 


36 BTU/hr/ft 2 /°F 


d. 


Thermal conductivity (k) 


6 BTU/hr/f t/°F 


e . 


Ambient temperature 


o 

o 

TJ 


f . 


Inside wall temperature 


1500 °F 


g. 


Heat generation rate (Q) 


432 BTU/hr/f t 3 



Theoretical solutions to this problem are readily 
obtained by methods outlined in Refs. 1 and 5. Comparisons 
of solutions obtained by finite element models and theoreti- 
cal means are given in Tables III and IV. 

A finite element model of this problem is constructed 
by assembling elements radially to form a segment of the 
cylinder. Since the problem is symmetric and no heat is 
conducted perpendicular to any radius, any convenient 
specification of element arc width, 0, and axial thickness, 

L, may be used. A small element arc is desirable since the 
isoparametric finite element used will provide a better 
approximation to a true circular arc. A large arc width 
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will introduce some distortion in the temperature distri- 
bution. Note that the sum of the forcing vector elements 
represents the total heat flux through the body in this case. 
The various finite element models used are given in Table II. 



TABLE II 

FINITE ELEMENT MODELS FOR RADIAL CYLINDER PROBLEM 









Model 






L-12 


L-2 4 


Q-6 


1 

l 

U 


Element 

Thickness (L,in) 


1 


1 


1.5 


1.5 


Element 
Width (0,°) 


2 


2 


6 


6 


Element 

Type 


Linear 


Linear 


Quadratic 


Cubic 


Number of 
Elements 


12 


24 


6 


4 


Number of 
Nodes 


52 


100 


80 


92 
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Figure 11. Radial Cylinder Problem Geometry. 



TABLE III 

RADIAL TEMPERATURE DISTRIBUTION IN A HOLLOW CYLINDER 



WITH 


THE INTERIOR 


. SURFACE 


MAINTAINED 


AT 1500 


° F 


Temperature (°F) 


Radius Theoretical 

(in) 


Model 

L-12 


Model 
L-2 4 


Model 

0-6 


Model 

C-4 


3 


1500.00 


1500.00 


1500.00 


1500.00 


1500.00 


4 


1188.54 


1188. 85 


1188.63 


1188.57 


1188.53 


5 


946.96 


947.34 


947.06 


947.01 


946.99 


6 


749.56 


749.94 


749.67 


749 .63 


749.62 


7 


582.67 


583.00 


582.77 


582 . 75 


582.74 


8 


438. 11 


438.38 


438.20 


438.19 


438.19 


9 


310.59 


310. 81 


310.67 


310.69 


310.68 


Heat Flux 


(q) 










Theoretical 


9.45 


9 . 45 


85.03 


85.03 


Calculated 


9.46 


9 . 46 


85.03 


85.02 
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TABLE IV 

RADIAL TEMPERATURE DISTRIBUTION IN A HOLLOW CYLINDER 
WITH INTERNAL HEAT GENERATION 







Temperature 


( °F) 




Radius Theoretical 

(in) 


Model 

L-12 


Model 

Q-6 


Model 

C-4 


3 


175.28 


175.34 


175.30 


175.30 


4 


173.00 


173.03 


173.02 


173.02 


5 


166.77 


166.79 


166.79 


166.79 


6 


157.13 


157.13 


157. 3 5 


157.14 


7 


144. 35 


144.34 


144.36 


144 . 36 


8 


128. 60 


128. 60 


128.62 


128.62 


9 


110.00 


109.99 


110.02 


110.02 


Heat Flux (q) 


Theoretical 




1.57 


14.14 


14.14 


Calculated 




1.57 


14.14 


14.14 
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3 . Chimney Problem 



Consider the tv/o dimensional heat transfer problem 
posed by a rectangular chimney with the inside surface 
maintained at 100°F and the outside surface maintained at 
0°F. The steady state temperature distribution may be 
obtained by constructing a finite element model consisting 
of the desired number and type of elements assembled to 
form a plate having the required geometry. The geometry of 

i 

such a model utilizing symmetry boundary conditions is shown 
in Figure 12 and the results of various finite element models 
are compared with a solution obtained by a relaxation tech- 
nique in Table V at the points shown in Figure 12. 
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Figure 12. Rectangular Chimney Problem Geometry. 



TABLE V 

TEMPERATURF DISTRIBUTION IN A RECTANGULAR CHIMNEY 



Point 


(x,y) 




Temperature (°F) 




Relaxation 

Method 


5 Cubic 
Element 
Mode 1 


10 Quadratic 
Element 
Model 


40 Linear 
Element 
Model 


a 


1/3 


45.5 




44.3 


45.2 


b 


1/2 


38.9 




37. 5 


38.3 


c 


1/1 


20.6 




19.4 


20.3 


a 


2,1 


39.2 


37.1 


37.9 


38.7 


e 


3,1 


47.7 




47.0 


47.5 


f 


4,1 


49.3 


48.7 


48.9 


49.1 
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c. FIRST ORDER TIME PROBLEMS 

Transient heat conduction in various bodies is a common 
form of the field equation that is easily modeled by the 
finite element process. Transient heat conduction in a 
hollow cylinder and a prolate spheroid will be considered 
here . 

Euler I, Trapezoidal, and Linear Time Shape Function time 
integration methods were employed. Each method yielded 
excellent results when used with an appropriate time step 
size. However, with the cylinder problem, the Euler I 
method was unstable for large time increments. The Trape- 
zoidal and Linear Time Shape Function methods produced 
excellent results in all cases. 

1 . Transient Temperature Distribution in a Hollow 

C ylinder of Infinite Length 

Consider the hollow cylinder of Steady State Test 2, 

Figure 11. Thermal diffusivity, a, is specified as 
2 

.25 ft /hr where a = k/pc. Initially the cylinder is at an 
ambient temperature of 70°F; suddenly the inside wall tem- 
perature is raised to 1500°F and maintained at this temper- 
ature. All other conditions are the same as for th_e Steady 
State Problem. 

An approximate analytical solution, which is valid 
until the temperature wave propagates to the outer wall of 
the cylinder, can be found by applying methods of solution 
discussed in Ref. 1. The result is: 
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The finite element models employed here are identi- 
cal to the models employed for the Steady State case. Tran- 
sient solutions obtained by the Trapezoidal and Time Shape 
Function methods and using the linear element model converged 
to the same results obtained in the Steady State Problem for 
corresponding model. Solutions with Quadratic or Cubic 
models and Euler I integration were not pursued to steady 
state . 

Tables VI through IX compare analytical and finite 
element solutions at various problem times for the 100 node, 
24 Linear element model used with the Steady State Problem. 
Comparable results were obtained with Quadratic and Cubic 
element models. The following abbreviations apply to the 



time integration 


methods : 


a . 


E-l 


Euler I 


b. 


TRAP 


Trapezoidal 


c . 


LTS 


Linear Time Shape Function 
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RADIAL TEMPERATURE DISTRIBUTION IN A CYLINDER AT .05 MINUTES 
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RADIAL TEMPERATURE DISTRIBUTION IN A CYLINDER AT .50 MINUTES 
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RADIAL TEMPERATURE DISTRIBUTION IN A CYLINDER AT 1.00 MINUTES 
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.00 70.00 70.00 70.00 70.00 70.00 70.00 70.00 



TABLE IX 



RADIAL TEMPERATURE DISTRIBUTION IN A 
CYLINDER AT 1.50 MINUTES 



Method 








Temperature (°F ) 








Exact 


E-l 


TRAP 


LTS 


TRAP 


LTS 


h (min) 






. 0005* 


. 001** 


. 001** 


. 01 


.01 


Radius (in) 














3.00 




1500.00 


1500 . 00 


1500.00 


1500.00 


1500.00 


1500.00 


3.25 




1240. 81 


1245.47 


1246.07 


1246.00 


1246.00 


1245.84 


3.50 




1009 . 18 


1016. 34 


1017.44 


1017. 33 


1017.32 


1017.04 


3.75 




806.91 


815.06 


816. 53 


816.38 


816. 37 


816. 01 


4.00 




634.79 


642. 85 


644.53 


644. 38 


644 . 35 


643. 98 


5.00 




220.69 


224 .17 


225.54 


225.50 


225. 39 


225.29 


6.00 




95.63 


95.97 


96.47 


96. 50 


96.42 


96.50 


7.00 




72.69 


72.57 


72.66 


72.68 


72.65 


72.70 


8.00 




70. 17 


70 .14 


70 .15 


70.15 


70.15 


70. 15 


9.00 




70.01 


70.01 


70.01 


70.01 


70.01 


70 .01 


*• 

.0005 to 
* * 

.001 to 


1 

1 


min, .001 thereafter 
min., .01 thereafter 









52 



2 . Transient Heat Conduction in a Prolate 
Spheroidal Solid 



An analytical solution for transient heat conduction 
in a prolate spheroidal solid is developed in Ref. 4. 

Results of a finite element solution to this problem are 
presented in Refs. 10 and 11. 

A non-dimensional finite element model using one 
cubic element or four Quadratic elements was constructed. 

The Linear Time Shape Function of time integration was 
employed giving results comparable to the results of Refs. 

10 and 11. A comparison of theoretical and finite element 
results is presented in Figure 13. 



D. SECOND ORDER TIME PROBLEMS 

The transient vibration of a perfectly elastic membrane 

is a two-dimensional problem that can be formulated as a 

finite element model of the field equation. 

Consider the square elastic membrane of length a and 

specified material properties T q , p, and v. T q is the 

initial tension of the membrane, p is the surface density, 

2 

v is the damping factor, and c = T Q /p. The partial dif- 
ferential equation describing the motion is: 



2 

c 




2 

c 






0 



where u(x,y,t) is the displacement of the membrane. 

An exact mathematical solution is derived in Ref. 2 and 
may be specialized to a simple closed form if an initial 
displacement of the form 
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Eh 



5 4 



CM rH 



at/L 



u(x,y,0) = d sin 



TTX 



a 



sin 



IX_ 

a 



and boundary conditions of simply supported edges are assumed. 

Two finite element models were constructed. One models 
the entire membrane as a 36 Linear element plate; the other 
as a 9 Quadratic element quarter plate imposing symmetry 
boundary conditions at the centerlines. Euler I and Finite 
Difference solution techniques were employed with a time 
step size of .001. Assume the following material properties: 



a . 


a = 


1 . 5 


b. 


2 

c = 


4.5 


c . 


v = 


3.0 


d. 


d = 


. 25 



Both time integration methods produced good results; how- 
ever, the Euler I integration method was unstable at large 
time increments. 

1 . Undamped Membrane 

The analytical solution is 



V 



d 



cos 



2TTCt 

a 



sin 



TTX 

a 



sin 



TTY 

a 



Results obtained with the finite element models are compared 
to the analytical result for the center node displacement in 
Table X. 

2 . Damped Membrane 

The analytical solution is 

y = d e (cos cot + — sm cot) sin — cos — * — 

co a a 
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TABLE X 



CENTER NODE DISPLACEMENT OF AN UNDAMPED MEMBRANE 









Displacement 






EULER I 


Method 


FINITE DIFFERENCE 
Method 


Time 


Exact 


Linear 

Model 


Quadratic 

Model 


Linear 

Model 


Quadratic 

Model 


.00 


. 2500 


.2500 


.2500 


.2500 


.2500 


.05 


.2378 


.2372 


.2370 


.2370 


.2368 


.10 


. 2023 


. 2007 


.2108 


. 2002 


. 2014 


. 15 


. 1469 


. 14 39 


. 1456 


.1435 


. 1453 


. 20 


.0773 


.0726 


.0763 


.0723 


.0759 


.25 


.0000 


-.0061 


- .0020 


-.0061 


-.0018 


. 30 


- .0773 


-.0844 


-.0787 


-.0839 


-.0784 


. 35 


-.1469 


-.1543 


-.1496 


-.1532 


- .1483 


. 40 


- .2023 


-.2090 


-.2041 


-.2073 


-.2027 


. 45 


-.2378 


-.2428 


-.2407 


-.2406 


-.2383 


. 50 


-.2500 


-.2523 


- . 2516 


-.2498 


-.2495 



where 














L 2 2 

0) = / A C 


2 

V 


and 


X 2 


= Tr/a 


Results 


obtained w 


ith 


the finite 


element models are compared 


to the 


analytical 


result for 


the 


center node displacement in 



Table XI. 
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TABLE XI 



CENTER NODE DISPLACEMENT OF A DAMPED MEMBRANE 







Displacement 






EULER I Method 


FINITE DIFFERENCE 






Method 


Time 


Exact 


Linear Quadratic 

Model Model 


Linear Quadratic 

Model Model 



.00 


. 2500 


. 2500 


.05 


.2389 


.2384 


. 10 


. 2105 


. 2091 


.15 


. 1716 


.1692 


.20 


. 1283 


.1249 


. 25 


.0854 


.0813 


.30 


.0464 


.0418 


. 35 


.0135 


.0089 


. 40 


-.0119 


-.0162 


.45 


-.0298 


-.0335 


. 50 


-.0405 


-.0434 


.55 


-.0450 




.60 


-.0445 




.65 


-.0404 




.70 


-.0340 




.75 


-.0263 





. 2500 


. 2500 


. 2500 


. 2382 


.2382 


.2380 


. 2100 


. 2089 


. 2097 


. 1705 


. 1691 


. 1704 


.1275 


.1250 


. 1274 


.0841 


.0815 


.0843 


.0454 


.0422 


.0457 


.0122 


. 0095 


.0128 


-.0128 


-.0156 


-.0123 


- .0309 


- .0328 


-.0301 


-.0411 


-.0427* * 


-.0406* 



-.0588 


-.0537 


-.0639 


-.0651 


-.0621 


- .0644 


-.0555 


-.0585 


-.0457 


-.0491 



* 

Changed step size to .01 

* * . . 
Changed step size to .01 at time = .5, Solution 

became unstable 
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V. 



CONCLUSIONS AND RECOMMENDATIONS 



The computer program presented in this thesis provides 
an accurate and reliable means for solving a variety of 
problems. The use of this program and efforts to increase 
its versatility are highly encouraged. 

The Linear Time Shape Function and Four Point Backward 
Finite Difference (Houbolt) time integration methods provide 
the most consistent and reliabje results. Investigation of 
other time integration schemes should be pursued. 

A. RECOMMENDED MODIFICATIONS TO THE PROGRAM 

A number of modifications to the present program are 
desirable in order to increase the versatility cf the program. 
Most of these modifications are straightforward and are the 
result of the simplifying assumptions made during develop- 
ment of the program. The recommended modifications should 
present no major programming problems. 

1 . Anisotropic Materials 

A facility for coordinate axis rotation will provide 
the program with the ability to accommodate problems in 
which principal material axes do not align. This may be 
accomplished by the addition of an appropriate subroutine 
to perform the axis rotation. This subroutine should be 
called from MERGE at the time of element formation. Minor 
changes to subroutine INPUT will be required. 
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2 . 



Time Dependent Boundary Conditions 



The ability to account for a time variation in the 
boundary value specified for {(f)} is highly desirable. In 
principle the modifications required are s traigh tf orward . 

The required modifications will only affect the solution 
subroutines TIME11, TIME13, TIME21, and T1HF23 if data input 
is accomplished by providing updated values of {<f>} and {<}>}, 
if required, in tabular form. The modifications will require 
attention to detail and a thorough understanding of the Irons 
technique as applied in the program. At present, boundary 
values of {<j>} are considered to bo zero; this effect on the 
solution matrix corrections and the effect of time dependence 
of the { F } vector must be carefully considered. 

Time dependence of material properties implies the 
recalculation of problem matrices at each station in time 
and possible iteration. Considerable computational effort 
is required and efficient programming techniques are 
necessary . 

3 . Non Linear Material Properties and Boundary Conditions 

A functional dependence of material properties and 

non linear boundary conditions such as the radiation boundary 
condition of the heat problem are real physical situations. 

The ability of the program to accommodate this situation is 
needed and feasible by various iteration schemes. 

4 . Problem Size 

The program is presently designed for problems of up 
to 110 equations with a maximum band width of 64. Practical 
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physical problems will require much more computer storage 
space. Modifications to accommodate more equations in core 
require a careful adjustment of the appropriate common and 
dimension statements and attention to the equivalence state- 
ments in the time integration subroutines. 

5 . Specialized Data Input and Output 

Hand preparation of input data is tedious and subject 
to numerous errors. For practical physical problems of even 
moderate size, automated data input is a necessity. 

Each class of physical problem will require data 
output in a specific form. Considerable attention should, 
be given to the development of various subroutines to 
tailor data output to a given class of problem. 

B. RECOMMENDED FUTURE STUDIES 

This thesis was concerned with the development of a 
computer program with the capability of solving a class of 
physical problems. This objective was accomplished and 
recommendations for improvements stated. Questions raised 
during the course of this development and beyond the scope 
of this work are: 

1. Convergence criteria of the various time integration 
methods, particularly the Euler I technique. 

2. Selection criteria for finite element type and, in 
particular, the type of element best suited to a given 
geometry or physical situation. 

3. Adaptability and usefulness of other finite element 
types . 
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APPENDIX A 



COMPUTER PROGRAM NOHENCLATURE 

A complete listing and description of all variables used 
in the program is not practical. The items listed in this 
appendix are common to several areas of the program and will 
assist the reader in a study of the program. For convenience 
items are listed by variable type.* The definition or func- 
tion and dimension, if applicable, of each item is given. 



A. INTEGER 


CONST/iNTS 




I GO 


Element type 




IPUN 


Punched output request 




IP ROB 


Problem type 




IPRT 


Number of time increment steps 
between print out of results 


to skip 


KBS 


Boundary condition surface 




KBT 


Boundary condition type 




KC 


Coordinate system of reference 




KINT 


Time integration type 




MRH 


Time step size index 




NBAND 


System band width 




NEL 


Number of elements 




NELBC 


Number of element boundary conditions imposed 


NGP 


Number of Gauss points 




NMAT 


Number of materials 




NNBC 


Number of nodal boundary values 


specified 
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NNODE 


Number of nodes 


NPEL 


Number of nodes per element 


NTSC 


Number of time step size increments specified 


NTP 


Number of times {<}>} is to be printed 



B. REAL CONSTANTS 



ABIGN 


Arbitrarily large number 


AKX 


[H] matrix property in the x direction 


AKY 


[H) matrix property in the y direction 


AKZ 


[H] matrix property in the z direction 


APZRO 


Calculated value of zero 


DT 


Time step size 


DTJ 


Determinant of the Jacobian 


FACT 


[C] or (G) matrix material property 


PHIREF 


Reference value of 4> 


TIM 


Initial problem time 


WFACT 


Products of Gauss weight factors 


WX , WY , WZ 


Gauss weight factors for specified direction 


XI , ETA, ZETA 


Local coordinates n / and £ 


X1,Y1,Z1 or 
XX , YY , ZZ 


Global coordinates x,y, and z 


C. VECTORS 
COL (32) 


Element level storage of the shape function 


ENDT (5) 


End time of the various time integration steps 


DTIM (5) 


Time increment for each time step size 


FE (32) 


Forcing vector at element level 


F (110) 


Problem forcing vector 
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IP (5) 


IPRT for each tine integration step size 


KBCOND (48) 


KBT for each element 


KEL (4 8) 


Element number for application of boundary 
conditions 


RHODE (1.10) 


Node number for specified 


KSURF (48) 


KBS for each element 


PHI (110) 


{(J)} vector 


PHIDOT (110) 


{<(>} vector 


PH I DDT (110) 


{<J>} vector 


PHIWRK (110) 


Working vector 


TIME (6) 


Time label for output 



D. MATRICES 



A J ( 3 , 3 ) 


[J] matrix 


AJI (3,3) 


[J] matrix 


BGC (110,64) 


[C] matrix 


BGG (110,64) 


[G] matrix 


BGM (110,64) 


[C] or [H] matrix, depending on problem type 


COARD (110, 3) 


Global coordinates 


CORD (32, 3) 


Global coordinates for element calculations 


DCOL (32,3) 


Derivatives of the shape functnon with respect 
to local coordinates. Column I is 3/3£; 

Column II is 3/3n; Column III is 3/3^. 


ENN (32,32) 


Work matrix for accumulation of [H] , [C] , 

or [G] e . 


ENW (32,32) 


Work matrix for element level calculations 


NCON (48,33) 


Connectivity matrix for each element, column 
1 is material number 


PROP (10 ,5) 


Material properties 


OUT (110 , 6 ) 


Equivalenced matrix for storage of output data 
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GP (6, 5) 
WF ( 6 , 5 ) 



Gauss points 

Gauss weighting factors 
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APPFNDIX B 



INSTRUCTIONS FOR PROBLEM DECK PRF PARTITION 

Although the procedure is straightforward, preparation 
of input dat£i for the program requires meticulous attention 
to detail. Errors are easy to make and difficult to locate. 

Input data preparation should, in general, follov; the 
steps outlined below before attempting to punch data cards. 
The use of the standard Fortran Eighty Column Coding Sheet 
is highly recommended. A sample problem deck follows this 
appendix. Integer constants must be right adjusted in the 
appropriate field. 

Preliminary preparation: 

1. Formulate the problem in terms of equation (1) and 
boundary conditions in terms of equations (2) and (3) . 
Identify the appropriate material properties. In the case 
of a one or two-dimensional problem, certain properties 
must be specified as 1.0; a specification of 0.0 will lead 
to a singular system of equations. 

2. Divide the continuum into the desired number and 
type of elements. Observe the restrictions of 48 elements, 
110 nodes, 64 bandwidth, 10 materials, and 48 element level 
boundary conditions . 

3. Assign node numbers as desired, observing the band 
width restriction. 
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4. List the element connectivity in accordance with 
the conventions of Figures 1 , 2, or 3. View £, n , and C 
as x, y, and z respectively. Any other convention will 
lead to serious error. 

5. List all nodal points and coordinates in any con- 
venient order. 

6. List the applicable boundary conditions and the 
appropriate codes and material properties. 

7. If applicable, select the solution subroutine 
desired and list the required time integration parameters. 

For the Houbolt method the specification of ENDT must allow 
enough time to generate the starting values required for 
the next step in time. 

Preliminary preparations are now complete; precise card 
punching instructions follow. Card format is given in paren- 
thesis followed by specific instructions where necessary. 
Recall that blank columns are read as zero. 

1. Title Card (10A8) 

Eighty columns on one card for any desired title 
information. Column 1 of this card will not be printed; 
if a 1 is punched, the output starts on a new page. 

2. Problem Parameters ( 1015 , 10X , F20 . 0 ) 

One card containing the following data: 

a. Cols 1-5 : Number of nodes (NNODE) 

b. Cols 6-10 : Number of elements (NEL) 

c. Cols 11-15 : Number of materials (NMAT) 

d. Cols 16-20 : Number of Gauss points (NGP) 
Specify 2 through 6; default to 3 
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e. Cols 21-25 : Element type 

1 = Linear element 

2 = Quadratic element 

3 = Cubic element 

f. Cols 2G-30 : Number of element level 
boundary conditions (NELBC) 

g. Cols 31-35 : Number of nodal boundary 
values specif ied • (NNBC) 

h. Cols 36-40 : Problem type (IPROB) 

1 = Steady state 

2 = First order in time 

3 = Second order in time 

4 = Second order in tire v/ithout [C] matrix 

i. Cols 41-45: Coordinate system of reference 

(KC) 

1 = Rectangular coordinates 

2 = Cylindrical coordinates 

3 = Spherical coordinates 

j. Cols 46-50: Punched output requested (IPUN) 

0 = No punched output desired 

1 = Punched output desired 

k. Cols 61-80: Reference value of d (PHIREF) 

3. Material Properties; One card per material for the 

steady state problem or two cards per material for the 
transient problem. 

a. Card I (I10,3F20.0) 

Cols 1-10: Material number 
Cols 11-30: k 

x 

Cols 31-50: k 

y 
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Cols 51-70: k 

z 

b. Card II (I10,2F20.0) 

Cols 1-10: Material number 
Cols 11-30: [G] matrix coefficient 

Cols 31-50: [C] matrix coefficient 

’ 4. Node cards; One set of cards for the steady state 

problem or two sets of cards for the transient problem. 

Each set has NNODE cards. 

a. Set I ( 1 10 , 3F20 . 0 ) 

Cols 1-10: Node number 

Cols 11-30: x coordinate or radius 

Cols 31-50: y coordinate or angle in degrees 
relative to the x axis (+CCW) 

(bis 51-70: z coordinate or angle in degrees 
relative to the z axis (+CW) 

b. Set II (I 10 , 2F20 . 0 ) 

Cols 1-10: Node number 

Cols 11-30: Initial value of 6. 

H r 

Cols 31-50: Initial value of <J>^ 

Leave blank if first order time problem. 

5. Node boundary values; NNBC cards required but omit 
if NNBC = 0 (110 , F20 . 0 ) 

Cols 1-10: Node number 

Cols 11-30: Specified value of <}k 

6. Connectivity data; NEL cards required for the Linear 
or Quadratic element and NEL x 2 cards required for the Cubic 
element (215,3213) 

a. Card I (215,2313) 
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Cols 1-5: Element number 



Cols 6-10: Material identification number 

Cols 11-79: Connectivity for nodes 1 through 23 

b. Card II (913): Use with Cubic element only 

Cols 1-27: Connectivity for cubic element 
nodes 24-32 

7. Boundary condition data; NELBC cards required but 

omit if NELBC = 0 (315,F20.0) 

Cols 1-5: Element number 

I 

Cols 6-10: Boundary condition type (KBT) 

1 = Volume integral of loading (Q) 

2 = Surface integral of loading (q) 

3 = Surface loss coefficient (a) 

Cols 11-15: Type of integral and surface 
applicable (KBS) 

0 = Volume integral; Use with KBT=1 only 

±1 = Surface integral over £,= ± 1 face 

±2 = Surface integral over n=± 1 face 

±3 = Surface integral over £=+ 1 face 

Note: KBS = 1,2, or 3 must be signed plus (+) 
or minus (-) to indicate the 
applicable integration surface 

Cols 16-35: Boundary condition property value 

8. Transient problem data; NTSC + 1 cards required but 
omit for the steady state problem. 

a. Card I, one card (2110, 10X,F20. 0) 

Cols 1-10: Number of changes in time step 
size (NTSC) 

Cols 11-20: Time integration method (KINT) 
Default to method 3 
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1 



EuJer I method 



2 = Trapezoidal rule method 

3 = Linear Time Shape Function or 

Backward Finite Difference method 

Cols 31-50: Initial problem time 
b. Set II, NTSC cards (I10,2F20.0) 

Cols 1-10: Number of time intervals to skip 
before printing output (IPPT) 

Cols 11-30: Time interval (DT) 

Cols 31-50: Total problem time to end of 
interval (ENDT) 

Data preparation is now complete. Cards should be 
arranged in order of ascending node numbers within each set 
or group and each set or group of cards stacked in the order 
prepared. Refer to the sample problem deck following this 
appendix . 

Several problems may be stacked and submitted if desired. 
In order to terminate a given computer run add one last 
data card with the word "STOP" punched in columns 1-4. 
Standard JCL cards and deck structure is shown in Appendix 
C. Print out of the program is suppressed by the "//FORT. 
SYSPRINT DD DUMMY" card; if a program listing is desired 
remove this card. The JCL card requesting additional output 
should be removed unless voluminous output is expected. 

Good Luck! 
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DO 800 I = 1 1 NP E L F0RQ280 

8Y = DC0L( 1.2 ) *DT J P0F.129D 

DO 800 J=1,NPCL FOR0300 

800 ENN( I , J )=BY*A<Y*DC0L( J , 2 ) + ENN( I , J ) F0R0310 



ooooo 
(Njcn^tin o 

OOOO^ 
aicccfu: cl 

oouuo 

LL U-LLLLLL 



OOOOO Oo^OO^OO^OOOO^OO 
r^c3 0'o*-«r\jrONj'Lnsor^coc> / ^^f\jro^Ln>or s - 
fo co o <- vt- <r <}- <t- vf <j- <r vr uom mir\ in lp« tom 
o o°o o o oo ooooo ooooo ^oo 

oc oc a cl oc rL cZ c< cl o' cL cL cL cl o' o' oC a' cL cLol 
OOOOOOUOOOUOOC'OOCJCOOO 

LLLLU U-U LLLLLLU-U IP LL tL U_ LL U . LL U_ LL LL LL 



O 

to 

Z 

o 



rvj 

O 



3 
O 
— I 
<t 

o 



H 




►— « 












K 






— 


H 




1 — 4 




to 


to 


O 




o 




CO 


CO 


UJ 




z 




< 


c 


Cl 




CD 




— to 


^ tO 


X 




o 


O 


to co 


to cD 


11 J 






1 — 


CD X 


CD X 


cc 


< — . 


>- 


o 


< - 


< 




o 


o: 




to o 


to “) 


X 


•» 


< 




CD h- 


CD h- 


h- 


♦— < 


o 


O 


x o 


X O 


►— « 


w 


z 


<1 


*> 


*■ •* 


j£ 


z 


3 


»- 


3 O 


3 O 




z 


O 


3 — 0 


UJ — cL 


UJ — OL 


to 


UJ 


CQ 


uj_jcl 


Cl 3 O 


030 


UJ 


+ 




CLUIO 


ZLUL) 


ZUJU 


> 




O 


zno 


-CL * 


- CL - 


►— « 


rn 


h- 


^z - 


<23 


< Z I 


H 


•> 




<r • _j 


h- -o 


b- *o 


< 


o 


LO 


h<cj 


UJ<CD 


UKO 


> 


*— 


z 


U 1^—0 


KJh" O 


rX 1 — o 




_j 


» — < 


NUJCJ 


•»UJ - 


► 3 * 


Ci: 


o 


<r 


I^J •• 


<NO 


< 5 


UJ 


o 


i— 


< *0 


h- *< 


b- -<x 


o 


o 


cl 


H<< — 


LU <£ ► ««-» 


LUC * 




o * 


UJ 


UJ 1 — 1 • — * 


► | 13 ~ 


►K 3 3 30 


UJ 


— if — — 1 Nj 


Cl 


►LUJLDw 


— .UJUJUI — 


^LUlJjLU LU — 


X 


LUOlU^Z 




►-t * ujcl _j 


X -n a 3 


x -o a a 3 


J- 


cl* cl< 


O 


x — . a_ z o 


— ZZO 


— — zz zo 




z~ z* 


Z 


— xz -o 


ox— *ld 


cx— •> o 


IL 


► ro * m 


»— < 




to— CP 


to — cDr^H- — <*- 


o 


•— 1 •> r— J CD 




Z UJ CO || “) 


ZUJ CD 11 O 


ZUJcO II O II < 




II »-< II || 


O 


n: cl o — . i— 


CLO-CD — H- 


Cl CL 0 »h^ -) || 


to 


►-< — O — 


_J 


O<o o 


□ <u o 


OCO — — 


H 


1 O 


_J 


lLX<rO || 


ILKO II 


ILKOmOO 


o 


OOO -z 


o 


tooO — 


z tooo— * 


Z t 0 -)O — o ►z 


3 


OUOmO: 


UL 


X O— 1 


CL >- -(•— < 


o:> r\j — i c\J h-*cz 


O 


000 — 3 




Ql — J — 1 *— < — 


OaT 33 Hw 


3 CL 3 3 r- 40 r -4 — 3 


o 


II 2 h 


UJ 


h J J _J 


1 — h- 3 3 3 hh Jj O Zl-O 


CL 


ONQ^UI 


X 


Z<<OOLUZ<<OOLUZ<<O II OZiuZ 


a 


Qcdqiuo: 


K 


luU(jooo:ijjuuoo^ l uuuo<oli'q:lu 




o 




o 


o 


o 




O 




o 


O 


o 




o 




O 


rH 


CM 








r— < 


r-H 


— 1 



Oo 



ooo 



107 



OOOOOr-iOOOOOO 


OOO 


OOOOOOOO 


ooo OOOOOOOO 


OO 




ooo 


— irxjo^j-tninvoxcooo^ 


(M rO <h 


invO^roO'O — in 


ro -tin vD X cOOOrHrXrO 


vt*in 




0 X co 


OOOOOOOOOO— irH 


rH , 




r— < <H *H » — < ' ^ (XI (XI (XI 


c\i M r\j c\J (\i r\i fM c 0 ro O 0 ^ 


roro 




rororo 


oooooonooooo 


ooo 


OOOOOOOO 


OOO OOOOOOOO 


00 




ooo 


CLClClQ-ClClClO-CL Cl Cl CL 


a, CL CL 


CL CL CL a- O- CL Q_ CL 


clclclclclClclClclclCl 


Cl CL 




D. Cl Cl 


iiiixirinxi 


XXX 


miiiii 


XXX XX XX XXX X 


XX 




XXX 


t/) 00 uo to 00 CO O) C/) oo UO </) 00 


uOL/)oO 


oooOL/n/)totococo 


l/) OO L/) CO tO (/) 00 00 00 00 CO 


00 L >0 




00 to to 


HOZ 








C —i 


o 


<X xo 


* r— 4 








“■"j — 4LU 




< 




h- M 


o o 


H* — X 't 


* > 








<t f—i ► 




H* 




LU # 


• o 


UJ — < O • 


— x- 












UJ 




M C 


CXI o 


M X OsO 


f— 4 — 








— nm 




M 


O 


•f 1 — 


t 


4- # OX 


M H 








rooc — 




4 


o 


o UJ 


# (\| 


O ~ 


x 












O 


o 


O M 


— X 


O X O -H 


< * 








roo^-H 




o 


• 


O + 


-H — O 


0 ^ —X 


K- • 








w Uh" 




o 


CO 


• O 


M HO 


• 0 •* 4c 


UJ X 








O ►O 




. 


X 


f— 4 O 


* MO 


' 4 O — *— * 


M -«• 








< — r- 


X 


f— 4 


— ■ 


— o 


< * • 


— O hX 


4 O 










» — 4 


■ — ' 


r-4 


■y- • 


►— c 4 - 


» M U 


O O 








— OO 


H- 




M 




LU h- \ 


— O' 


O O 








<M 1 « 




— 






M UJ — 


-H 4. <X 


O • 








mo^ 






c 


> * 


4* M-H 


O h- 4L 


• 0 > 








r— * 


O 


> 


y~ 


* — 


O 4 » M 


* O LU *— 4 


4 








_j — 4|— 






UJ 


c •— * 


O O^; 


<HO mq 


— O 








O — Q 


CO o 


< 


M 


H* X 


0 r>< 


h-O . 4 - CO 


* O 








USD 


< O 


h~ 


4 


UJ -K- 


• Oh 


UJxt- H O • 


— O 










► ro 


IU 


O 


+ < 


-H • UJ 


4- •*— O r *H 


H * 








— COX 


O ro 


4 


o 


O v— 


HM 


Ov0 4c O 4- 


> —4 








ro ^cl 


a: ro 


o 


o 


O UJ 


* — 4* 


Om— •*— 


* — 










m ro 


o 


« 


o + 


— 41 O 


O — — « hX 


< # 








(\un— 


Cl — ro 


o 


•— e 


• o 


-H — O 


'•hX — X 


y~ — 








ro •'O 


<r r\i ro 


• 


— 


r— 4 O 


>- —‘O 


hO^ r — 4 


u » *-» 








- 0-1 


-— ro 


r-4 




— o 


* > • 


— —O 


4 X 








i . — * . — i 


u ,0 ro 




— 


* • 


< -H 


u • X -hCO 


0-4* 








O — — 


aicC ro 


# 


^-4 


O •— J 


L— C — 


— • 1 >- • 


00 — - 








OCL^ 


c^h- ro 


— 


>- 


-HO — 


UJ H «■ 


-H | O 4» -H 


otx 








ooa: 


*— * .-xxxrO 


r -4 




Xo^- 


4* UJ — 


x— n < + 


• * 1 








•-CCZZ 


X — OOorO 


X 


< 


•W* • — 


O 4* *— * 


■JJ- <X O — H ^ — - 4 


^Oo 








'-'CL*-' 


cl stnOOrO 


•K 


k- 


*—• CO *— « 


00 OX 


• O UJ X 


— xO 








— <M -X 


r-ooon 


— * 


UJ 


XmX 


OOO^ 


XUJH 4 - 4- 4i 


•X —0 








-J ro— a. 


k— m • • .ro 


X 


+ 


+ — * 


• OO — 


4* M . .0—4 


O H • 








UJ +CQ * 


oO^^^ro 


4 


O 


Oo •— 


H 1 • X 


o > -«■ OOO 


OX— 4 








Cl f\l vt'-' 


<XOC - | | • 


o 


a 


oox 


— CO ^H 1 


OCOMOO 


Ok— 








~Z. ro— o 


LLOO •* •* I 


o 


o 


Oo 1 


* x — 0 


Oh- 0 — • • 


. 0 > 1 - 








► — co^ 


— *0000 * 


o 


• 


. .o 


-4—* O 


• LUOhhh 


C'OO 








< ^ Cl— (stM ►oooo 


. 


r-< 


-H(\IO 


x 0 *— *o 


hN • M — | 


II 00 








Zo- 


o • ooo 


r— 4 




- 1 O 


11 0 x • 


— 4- O'* II — 


— *o 








lu u j cc h- 


*■ <x st —4 9 • ro 


— 


* 


II r4 • 


*— * O | * — 4 


II < II < — < 


»— 4 Qx • 








M — »-D-QO *— | 


II 




— M -H 


-H • II — 


— > — • — ^ 1 — • — « h— 


m + cr* 








•^M-^ O »— < >- >-. * | •* ro 




X 


-H^f ^ 


M — 4 — * 


HUJhIUNUJ 


* •—* 11 




00 




< i 


h^oo -oco 


f — 4 


II 


MC II 


► | — H r*4 


M-tf MM ►M 


— ‘X — 








t-O »-f\JX 


^ < a" ooo co 


M 


— 


•“H — 


hhM> 


► < - 4 


>-3 — 1 — f 




0 




uj ^cxicoclo •‘noooo'i 


•- 


f— < 


— H LU — H 


>- M ^ II 


HhHO>-< 


► — M 0 




►— 4 




•‘Xro— - 


0X0 • O .ro 


f— < 


M 


>-MM 


+ ¥ -J— 


>uj>o • j— 


— 'X •“ nT 


O 


K- 




— • 1 '-'LU — 


COX - — » *—4(0 


>* 


»■ 




-H C >~-H 


• 4- ^O— ‘LU 


X-if— 4 . 


O 


O 




X< OLLO 


+ <x — | —4 i ro 




f— 4 


-H ^ f— < 


X \ — ^ M 


r—i * — 1 *~H • X M 


0 


►— 


2 




— ' — OX + r-K 


— * ► *ro 


f— « < 


X 


x>> 


“UJC - 


xxx — 1 - 4- 


<X * X 


»— 4 ► 


O 




uj O^^sfo^oooro 


X 




►-U- ^ 


<MhH 


•* •“ — <<! 


K- ^ — 4— * 


4 — 


LL 




clcoo<m— 


^H-xoOoro 


r 


r— 4 


<c<r 


H* 4- UJ> 


<— «<X4L h- h~ 


UJ— iX —4 


(MO 






<■* ro*— 


•"Oooooro 


< 


X 




UJ » — 4 M ► 


L — X h- —• U 1 UJ 


mO » r^j 


—4 ro 


h- 




I-JS •lO o' • • .m 


H* 


•- 


UJ UJ UJ 


MX ►C 


LU'—LU-HM’X 




x - 




(\j <r 


to < x rsj cl — i x o — * —4 — i co 


LU 


< 


M 4- M 


- <h 


N| M> . C 


< • 1 — <1 


— O 


UJ 




LuaTro - 


— icloMxxco 


M 


h~ 


4 ► 


<<hiu 


^O c, h- 


1 — (\|UJ h- 


CO (M 


S 


W^4 9—t 




— -J X>~M • 




UJ 


<x< 


|— LU fvj 


< O <T <T 1 — LU 


LU | MUJ 


1 ^ 


UJ 


II II 


2 5T30ILlZn00\ 


COM 




UJuj - » 


h- O k— h- m 4 


► — < •* M 


-JO 


-J 


MX 


^hNZ^covOorcXcLO 


1 — o * 


LU *— « LU 




LU .UJUJ ►*— « 


— X^- 4 


LU —4 


LU 


• — « * — * 


UJr^CO -OOO^ 


UJOC 


*■ X - 


0 — -HXX 


-QX ► 4 — X 


xo 


CL — 






002 *— 


►Z^OOUh 


•* 


• 1 — 




oxx*— ^ 


M *^ M OX^ 


v^o—O 




CC 


OOO 


O'-'O— IL 


-oz 




CO LU 


X-tt xo — -u x> 


x-^ xo — — 


XO>o 


wO 


< 


II 00 


a:-J5:ro - 


ia?:uj< <<< 


Xx — 




• 0 •— X M 


— — —oox 


5*OS • 


II 1 — 


LU 


) r— 4 r— 4 


cocls: 


*— * X 1 — f — ) — » | — 


w< 


— -J 


UHt^OXOO 


UHV «0 — 


0 «o— 1 


0 




O 


-j^-OroroH’CJ — <<x<<r 


-J — <LL 


CNO\LL-ft LULL 


UMUhu. 4t 


LLfsJLU — 


00 




OOO 


</)*~*o — *• 


QUQQQqQ 


LL MO 


U-Jt LL- 


— 0000 


LL LL — Q *-H *— 1 


10 104 


—0 


-J 


m OO 


— KMCO 








r- < 


H r-H ^H(\ 


1 *—4 —4 






0 



o 



u 



u 



u 



ooo 



OOOOOOOOO 

ChOrHfMrOvfin^r- 

fOst 

OOOOOOOOO 

O-CLaaaaaaa 

iiiiiiiii 

</)</) oO LO l/KO LO oO oO 



ooooooooo oooo 

cOO'O-- fCMrn^invOf'-coO'O 
vtsTininiALninim^^ run 
OOOOOOOOOOOOO 
CL CL a. n. CL CL CL D . CL Cl D a CL 

X X XX X XX X X X X X X 

oOlooooooololOoolo ooojoo oo 



OOOOOOOOOOOOOC^O 
r-4 rg o ^ in o r- co o o nj rv> 4* m 
v 0 vO'O'ONO'O>O'Ov 0 r^NNh-r'-N 
OOOOO^ OOOOOOOOO 
a o a a a a a a a a. a a a <x a 
iiiimriixxiii 

OOLOLOoOLOUOoOl/)OOoOl/>OOLOOOoO 



r— 4 fH , — 4 

XIX > 











► ► ► 
















r-4 —4 r-4 






r— 1 


^ f— ( ^ 








^XNX 






X 


*—< Xj *— < 








Xl - 






► 


M -X— l 




• — < 




** r— 4 *“H r— 4 






^ f— 4 


► H ►>- 




*- 




r4X>* Xl 






- — . *—4 — r^j 


•* 




x 




X - •* •“ 






r-4Xj r— 4 ► 


> -xl-4 




t— 4 










X 1 X 


•**— * h X 




► 




-4hXh 






♦ •>+— 


hXh - 


00 


in 




XUJ *-iu 






.— < *— 4 r— 4 LU 


X 


7Z 


r— 4 




►xl <1 ► 


COO cm 




N>XN 


•■<1 *■ xj 


O 


► 




< 


•—4 CxJ ' — 4 




► r r » 


<>-*-< - 


►— « 


cO 




h- <IU X 


» *" » 




< — << 


h-UJXC 


h- 


r— 4 




LUf—Xl •• 


sO • 1 




hXhh 


UJxl -h- 


O 


*• 




NUj #- < 


r-l *— ( r-H 




UJ ►UJUJ 


XI *-<tUJ 


2! 


x 




► ►< h* 


► *• t- 




XI r •» 


-<th- - 


o 


*• 




•X » i | UJ 


vOCOO 




♦ i — • — * • — * 




X 


in 


CM 


hXiUrsJ 


rH 




<uj >: x 


hUJNlX 




•> 


•* 


UJw^-w 


f» ^ r 


vT 


1 — xl * ^ 


Liwwv^ 


H 


CO 


H 


•*000 


CvjNf o 


► 


UJ X 


- 1 1 — > 




•* 


» — < 


HOOO 




h— 


XKX. 


— LLLLlL 


ID 


r— 4 


•- 


XllU-U. 


• • 




— I— UJO 


xooo 


s: 


W 


*— 1 4 ►*>< 


—oon 


to 


CM t» 


X LU CM u. 


w II II H 


UJ 




cm 4- •— 


O II II II 


UJ 


c\l 4* *— 


— wwQ 






00 


* — II 


^ o — * ^ 


o 


r r-«* •— < 


— 11 


— <X XrvJLL^cNJO 


UJ 


UJ 


p-— 4 r— 4 _J 


X X xi u_ * — f r\i ro 


o 


Hr- 1 H 


zr^ooo- 


+ m mm II *-*-*- 




o 


II | o 


►— <»— «*— 1 n * ► 


2 : 


II 1 CM 


mmmugLLH 


-J- — — 


o 


o 


NlNl o 


_J_J _J 




rMiM OCM^H || || || 


OX^MJOOOO 




3 


i »— i 


>*XrvJJOOU 


UJ 


» — « >— < > — « 


4- +OOX>nJ^^^(M 


UOQQQuUUvJ- 


K 




~0 


^OOOOOOO 


o 


WsJ- 


(MG' | OOQrjg-O'U 




< 


a 


o-?:- -f o + o^o^ctu^ M *-' 


*— • 


o* + 0000 ^ 0 :^ 0 : 000 ^ 


II □□OhhwwwG 


(X 


UJ 


OCM ' OxOOO^ — —w 


00 


co in cm*— m 


• — * > — < O CO O O * — ‘ 


JOuU^JJJh 


o 




C\l •— 4 ' O CsJ 


^UUL) w ) — 1 — I 




|| f\J r -4 »— « CM 


II II •-•OOOwww ^ ) 


O II II II JOOO 


< 


a 


II II II 


II II II II JOOO 


o 


O' 11 11 


11 II 11 11 


UHrg^OljOtJO 


3 


o 


OMhXQXHHr4GUUu 


*— « 


OO^HOUUZhh^OOOu 


h X>'N|(JOQQO 


C3 


o 


OmHMQ 


*— X XMOOQO 




mQmhQ 


mmmX>NuOUQ 



o 

o 



o 

CM 



O 

O 

CM 



UOOOOO 



ooooo 



109 



oo n oooooo 

h-r-h-r-cococococo 
oo ooooo^o 
a cl cl cl cl cl cl a cl 

iixmiii 

bOOOOOoOLOOOuOOlX) 



OOOOOOOOOOOOO 

Lnor^coO'o^fMrnvtinsO^ 

COCnro O 3cO0'G-0'0'C' r J'O'n^ 



CLCLCLCLQ.CLCXCLCLa-O.CLCL 

IIIIIIIIlIIII 

OOOoOOOOOUOOOLOOQOQuOLOLQ 



000 ^ 00000 oooooo 
coc^o^'xi-n^iAvor^coO'o^fxi 
^Q'OnoOOOOOO°-<-<^ 

O O r— 4 ' — I •—■4 r- < <— 4 r— 4 r— 4 r-— 4 r-4 »~4 , — 4 ■— 4 

CLO.CLQ. Q-CLCLCLCLa_CLaCLCLCL 
III1IIIIIIIIIII 
L0 l/)L 0L0 o 0CA;U0O0C/>O0O0L0 o 0L0 u 0 



-> >- r-^f — . — » — * — 1 

•* I — (>-r-lX 
>-*“-1^-4 **X **X •* 

^rs|rvj< I -i ►< 
i— l •* I I — ► X'—U- 

fslr- 1 ^lUr-4 ^>-UJ 
* >- <3 X >- <t •* 

<t *■ I — + ► I — ►— * * — • 
H<LUh-.^LUXX 
LU t — MX X - •* - 
fslLU *. ► *<<<J 

» — « • — < I — h- V— llj LU UJ 
XXUJUJLUNlfsJN 

2 : 

0000000 o 

LL LL LL LL LL IL U 11 

OOOOOQOO 
II II II II II II II II 

(Mrn^HC\irQrHf\jra 

cm cm >4* O' O' O' o 

UUUUOUUUvf 



U1 

o 



o 

2 : 

O 

LL 



7" 

UJ 

UJ 

-J 

UJ 



ww— O 

-J-J-J-J_J.J_J.-Jh- ~ 

00000000 CO 

OOOOOOOOO o 
QOOOQQDOo o 



o 

in 

CM 



o 

00 CM 
CM— I •* 







— 


►m 










r — 4 r — 1 *~4 


O -r -4 










XX>- 


CM vO - 










»■!'*> 


-CM 00 










«— 4 r — 4 *— 4 


CO ^r - 1 






0 




— «>-xix 


r -4 •* 






0 




xl 


► r — 1 I s - 






9 




•**—4 r -4 r — 4 


m * H 






h- 




-HXXX 


O 






CM 






X 










•‘<~ < 


XxM 










HhXh 


XX 






CM 




XUJ •‘UJ 


► 






•* 




•*X <t ► 


O' CM O 






H 




< ► X •— • 


CM CO ^ 






CM 




|— <IUX 


► •* 






X 




UJ | — Xl *- 


cm in in 






O 




XILU *.< 


CM CM 






r>H 




► •*-< h- 


»»■ *■ 






•* 




<M-,^UJ 


O' CM - 4 " 




O 


X 




HXlUXI 






— 


► 




1 1 1 «— » >->■» 


p* + + 




► 




h- 


•^OOO 


CM m CO 




h- 


•* 


►_ « 


•—OOO 


»~4 




►—4 


r — 4 


► 


XLLILIL 


• • 




+■ 




r — 4 M 


— 00O 


LO 




XXI/) 




CM + ~ 


0 II II II 


UJ 


CM 


x * — * *— • » — * — ■ 


l/) 


•^ II 


^ 0 *- 4 * *■"'* r-'* 


O 


•- 


| II — X 


UJ 


r ~4 r — 4 _J 


X X X LL r-( CM cn 


O 


r — 4 


r -4 0 ^O_J XXX — 


0 


II I O 


» — « ►— 4 >— 4 ll ► » * 


2 


II 


II -<—10 


0 


XI xl 0 


— — — 1 — 1 — 1 




XX 


*— • 1 1 O*- 1 ^ — — — 0 




►— « 1 — * » — 1 


X>xf JOOO 


UJ 


» * » — * 


1 — i • — < ►— < » — 1 + -f X >X!cl 




— O' 


*— 1OOOO0OO 


0 


* 


►— * 4—4 — i^DODl — 


cl 


0-* + o + cl cl cl o 




cmocm 


O-ft V Or-iO + cCcCcC II 


UJ 


00— 0 X 000 —' — — — 


00 


r -4 in | 


inon^M^xOooo 




01 cm O rO 


— 000 — -J—l-J 




11 cno'OcncMCMcn h ►-« — ooocl 


Cl 


11 11 11 


II II II II 0000 


0 


cn 11 


11 11 11 rn 11 11 11 11 11 •— 


O 


O^hXOX-iHr-iOuOO 




rioOXOi/)HO^u'X'-i^^x 


O 


QmmmQmXXMOQOO 


2 : 


— O — 






O 


0 










CO 


0 












0 









000000 



00000 



110 



OOOOOOOOOOOOOO 

cO vt ir\ sO m rr> C' o —< oj co ^ uo o 
r-t ^ ^ h ,-< r-i ^ rsj cm rsj rsj *\J <\J (M 



<xa.a.cLCLQ.Q.a.Q.a.a.a.CLa. 

xxrxxxxxxxxxxx 

oo oO co oo oO co oo oo co oo co co co co 



w— -4 ^ 




X -» 


r— 4 


•“ M 


>- 


* -">■ I 


«««, * 


M «— 1 i - 




•>-^>^ 


~>-^X 



>.r^^Q *“ M *»>-r— 4 **X * 
•‘M >- Cl —4 ►r-< •'rvjr-l eO 
r-4 | •> ♦— < M * — < M Q | X »""* OH 

M ^ r-H X »•>• | CC ► — 

rHXf-C - “HHO -x 

0 >- - ^OQI>-n:Oh 

cl ► o < ►-< r£ »— -+-+c<: •* 

hQq>I h m •‘OX — < 
Xc^— 'iiJ^XX — c^^-Xh- 
•- »— h- X — •‘K- UJ 
-Ih “<I * •* -X< - •* 
< 1 — •«<}— — •— 
h- * < h- U ) h- H- h- •‘LU X X 
tU»— UirvJlUujUJ*— •*•■•* 

MX LU •* *• *-MMX<<< 



<< — XXX<<<UJIUUJ 
h— X — l — I — MMM 

LLltU **X>- >UUUJtU — — — 

*M <1 >Z7£ls -' — ^ ^ X 

•HUOO>-X>-^X>- 

x<luu, u.u-> ^rooo 

wj—MOOOt )UUu. U-U 
xiu— ii ii ii iulu-DOO 

'-ocQ ii ii ii 

ILM O r-lCNJCO II II II ~ — — 

II t)U- - - ^^(vjro 

— ul ii _j _j .-J ^ og co - » •* 
Jll-'OOO •* “ -rococo 

O^TNr— — ' — « *— » — ►— I 

• — •• — . • — • w w w w- ^ ^ ar 

^W^JJ-JJJJJJJD 

JJJOOOOOOOOOHQ 

OoOUOUoOUOOOUJ^ 

uauoQOooooQoa:Lij 



oo 

iTwt 

CO 



ill 



SUBROUTINE JACOB ( NPE L , A J , DCOL , CORD , A J I t DT J ) JACOOIO 









OO 




00000 




003 




OOOOOOOOOO 




0 








r\jro 




vt insO r-00 




O'O — 1 




(MOvtir.sOr-cDO'O— « 




CM 








OO 




000 00 




O *— 4 1 — * 




r ~4 r—i f— H f H f-J r— 4 r-H f\| C\I 




CM 








OO 




00000 




OOO 




OOOOOOOOOO 




O 








OO 




00000 




000 




OOOOOOOOOO 




O 








<< 




<r < < < < 




< < < 




<<<<<<<<<< 




< 








Oo 




0 0— > 00 




030 




000300 3-)0 0 




O 


y- 
























2 


H- 






















< 


O 












CMfO 












UJ 












< V 










•-H 


CL 












3 CO 






> 




21 


00 


















3 




O C 


LU 












0 






2 




IJJ 


a: 












— < 






O 




J — 














CO 1 










LU 


X 


















CO 




Q 


y- 












r-H *— < 






2 






»— « 












*— «■ *» 






O 




00 


3 : 












0 3 






y—i 




H 














<— ' 






y- 




»— i 


00 












+ 3 






< 






UOLU 












— < 






3 




LU 


LU > 


















3 




h- 


y— *-« 


















O 




3 


<h 




— * 








COCO 






3 




Cl 


2 < 




CO 








— » 




OOOO-JOOOO 


< 




>7 


.-> 




•> 








O CM 




l~ ►— ft — » — ft — H-H- 


U 




O 


O — 




CM 








< — 




OOOOOQOOO 






O 


cr or 




CO 








* 3 




SMVMSMSSS 


2 






OUJ 












— < 




^ 


O 




Q 


2 OO 




— J 








COtt 






> — « 






< O 


K— 


O 












CO CO CM c O CO CM CO ro CM 


H- 




< 


*-< 3 


2 


O 








<MCM 






►—4 






cnp^u< 


IU 


O 












CM OJ C\J rH w—4 »~H r-j r— < 


O 




2: 


o<<~ 


00X 


*- 








0 CO 




WWWWWWW^W^ 


2 


— 


< 


O ►—•CDL— UJLU 


— <. 








< — 




OOOOOO OOO 


a 


LO 


» — < 


<cQOa'hJ 


CO 








3 




<<<<<<<<< 


0 




CO 


oc>3<<iu 


*- 








— < 




± 4 k % K * i * K 




r* 


O 


uoa 


2 


CM 








CM | 






>- 


O 


O 


xiu< 




CO 




•» 




► — . 




CnJ^-<^hCM^-(— <C\J r-(*— c 


Cl 


y— 


< 


— ‘XOLLLLOUJ 










» — * CO 






< 


0 


— > 


c^h oocc a_ 


a 








W » 




co ro ro 10 co cm cm 'M 


a 


** 




h- UJ 


0 


or 




0 




0 3 






2 


0 


LU 


<U_XXXOoO 


0 




a: 




<W. 




000000000 


3 


CL 


X 




OUJ 


0 




0 




+ 0 




< < < < < < < < < 


a 


O 


1 — 


CiToT 


O 


•- 




0 








1 1 1 1 1 1 1 t 1 


CD 


LJ 




2UJLLH- h- 


OO 






* 




CO A 










zr 


< 00 O<<<T 2 T 


-SCO 




— 




► — * 




P^COCMCOfOCMfOCDCM 


O 


3 


ct: 


— Cr: 2 X 0 


rv} 




* — « 




cO *M 






V — 


O 


0 


CO LU V— 


on 


1 CO 






V- 


W * 




CO p"i P^ cO 3 r\j cm (M 




C 3 


UL 


o>zroro->n 


a— 






2 


O CM 




WWW 


uo 


O 




oz< 




*- > — « 




— 


< 


<W. 




OOOOOO ooo 


2 




0 


<mZXXQQT 


X-> 




_J 


2 


* 3 




<<<<Z<<<<< 




O 


— 1 


*0 — 


H-LU 


1 < 




0 


*— 1 


— < 


UJ 


<***#*<*< 


< 


< 


hH 


mzjj 


CD 


< •* 




0 


X 


CM fc 


oO 




► - 






ro a: lulu 


X 




2 


0 


Cl 




07 


(M^r-4CM^|rH(\J^-«^-< 


o: 


3 




xujclcl 


3 


CO 


< 


+ 


UJ 


CM 3 


UJ 




UJ 


UJ 


IU 


X 1-22 


2 


CO * 


» — • 


-J-' 


y- 


w , — . 


> 


CM CM C \J , — 1 <— * —* r— 1 r— < 


a 


CL 


2 : 


coiu 




# CO 


CD 


LUO 


UJ 


3 COM 


2 


w'-'wWWwwWw 




2 


M 


m 0 11 11 


II 


3 — 


O 


OCL * 


O 


< w ► 


►— 4 


OOOOOO ooo 


0 




H • 


II 11 0 3 


O 


< *0 


O 


comox^ 




* 3—1 




<<<<<<<<< 


2 


CD 


3uj 11 ►-•oclo 


UJ 


UJ< 


< 


r » O * 


U) 




UJ 


« — Wwww — .WWW 




CD 


Di/)“ 5 “)hOiJ 


X 


QZ 


”0 


1 — 1 r— < • r— » O 


X 


3 | O 


X 


till 


3 


O 


nroz<<ooo 


2 


2 




II ll O || < 


H- 


- — < 


1 — 


II II II II II II II II II 


O 


O 


CD UJ 




>-o 


UJ 


►~»0 II II 




3 CM* 






u 


< 


3 > 






• »— 1 


X 




UJ 




UJ 


r-A *-4 CM CM^M p^i CO CO 


3 


O 


O02 






Ooo 


h- 


00000 


y- 


OC03 


h- 




0 




> — < 


1 




— 2 




00 **0 *“ 


3 


<w * 


3 


r-1 CM CO r-l (M CO CM CO 07- 


LL 


> 


00 






OUJ 


z~ 


LA UO 1 LTN •—* 


Ol 


II O CM 


CL 






or 


~CS 


1 




Ci_X 


or 


w w 


jr 


o<-* 


2 


» — < ►— < ►— « < ►— « » — « ►— * » — < > — « f— 


UJ 


H- 


12 






X~ 


0 


00003 


O 


h* 3 


O 


OOOOOO OOOLU 


X 


2 


H< 






— Q 


LL 


oo<o< 


O 


0 — < 


O 


<<<<<< <<<c^ 


h- 


UJ 
















»— 1 CM 




















0 
























0 
























in 















ooo uou ooo uoo o 



112 



FORM THE JAC03IAN 



OOOOO 

ro<fir\vOc— 

(N](\J<NJ(\JCNJ 

ooooo 

OOOOO 
<<< < < 

OOOOO 



OOOOOOOO 

COOO^NrOvfin 

(NjNrnn^rnrnm 

OOOOOOOO 

OOOOOOOO 

<<<<<<<< 

OOOOOOOO 



o 




(NJ 


ro 


(NJ 


*“ 




•- 




•> 






ro 


r— 4 


«— 1 




t; 


W 




w 


o 


7! 


o 


o 


o 


cc 


< 


< 


< 


< 


o 




•# 


* 


* 


o 


►— < 








* 


V 


m 


r~4 


r-4 




a: 


► 


9- 


*■ 


* — i 


LU 


<NJ 


ro 


(NJ 


•> 


h- 


W. 


w 


* - 




UJ 


00 O 


o 


o 


w 


o 


*:< 

i 


< 

i 


< 

1 


o 


UJ 


* i 


i 


1 


o 


o 


Oro 


CO 


(NJ 


o 


< 


O - 




9 


+ 


ll 


Oro 


CO 


(NJ 




or 


► W 


W 




UJ o 


o 


oo 


o 


O 


OC3L - 


oo 


o< 


< 


< 


cocoox.-* 




CO* 


* 


* 


*- t-O •* w 


UJ 


► 






r~i —* »HO 


X 


O (NJ O •— 4 O < 


II II O II < 


h- 


o -a 


l rO 


1 9 


HHO II ^ II 




h- C\J O *— < O •— < 




UJ 


r — ( 


Iwr^ 


Iw 


0 O 000 


H* 


O 


o 


OX 


oo - 


o 


0<G<0<Q^ 


O \0 *— < v£) »— • 


Q- 


H- || H 


tt h- 


II o 


W w 


51 


o 


o 


Ol — I 


OOOOo 


O 


OhOHOhLU 


oo<o< 


o 


OQOQCQchi 


O 




O 


o 


OO 


O 




o 


O 


OO 


o 




I s * 


00 


a^o 



oo ooo 



113 



o 




0^00^00^0000000^00^000 


r— 1 




r-iCNJOD 


LnoNanoo-irvjfDvtinOhccG'O-i 


o 




OOOOOOOOO*-^ 


> — < ' — < r—{ w — 4 r— < *— < *— < r~4 fN J f \J 


o 




OOOOOOOOOOOOOOOOOOOOO 


OO 




OO OO OO oo LO 'XJ OO OO !XKO LO '^0 OOlO OO OO oo OO L/; OO 


o 




qdqqoodddcodqooodooqd 


-J — 




— 1 -J -J — 1 -J -J _J _J — 1 _J -J — 1 -J -J -J _J -J 


-J-J-J 


-J 


h- 

* ~ -J 




* 








* 


* h-cn — 




* 








* 


* or-J<r 21 — 




* 








* or 


* <— — 1 moOX 




* 








* LU 


* UJ— ID LU*~‘ 




* 








* CD 


^•200 dcSCS 




* 








* 2 


* — 2 LU Y- 




* 








* < UJ 


*-J~<r 2lU< 




* 








* lu >r 


* -M-2orx2 




* 








* >-X oo 


* LU-O^DH- 




CM* 








* CDh- K- O 


* x uJDh a: • 




X* 








* 22 


* HHQ' JiuO<^ 




o* 








*oor iu<n- 


* O Oft'Z-JH 












* UJ LU »-h x *-< 


*u^3<0 




4f 








* 3ih- o x 


* GO oO OLU 












* OLL UJ 


•If O 2 1~ h- 2 2 2 




<* 








* -J<H-U.X 


* * or oo —J 2 < uj 




5:* 








* -J — ILL O 2 


* < Cl oT ID ID — -J 




* 








* O ►Oiu^O 


* O — uo _j crT uj 




*- * 








* LL — -JOOTCL 


* XLUUJX UJOH 




-J* 








* CD O ID 


^ M Ia' oro -J 




o* 








* tu or 


■K QCJ-CJW H < 


51 


o* 








*2 2LL O CiT 


•if h- \— X u If— »— • 2 


UJ 


X* 








* a ii — o a 


* <O^H Xoo2CD 


h- 


o* 








* — • o' h- 


* x \ — i-ariDo 




OO* 








* — uxoo 


* 2UJ2 ~ < 


— XX 


* 








* DXO^hUJ 


* f-^aQ^OU-tu — 


or ooll 


UJ* 








* I- — OTclO> 


*2 2 X O 


h- h- 


h- * 








* ‘ CO | — UJ 


* UJ *h-oo<T LUh- 


< a i oo O 


<* 








* H-* ID OOO 


* *-<^oOh- X h- 


2-XXCiT 


=)* 


o 






* oO < t/; r-_ »— < 


* UcODZtOKU.^ 


► -O0UJ 


o* 


o 






* CD— OX 


* •-«— 2 x 00 ox 


I- M 


<1* 


o 






* ID UJO<K 


*u. 2 iu2 3 


Z2 liJ 


oc* 


• 


o 




* ^XlUJO 


* X II K-Ujo— ‘00 


UJ^XU- 


o* 


o 


o 




* LUI— O0-J2 


* uj -h_jo i~o 


— H*0 


I-* 


II 


oo 




* Of- O ~ 


4c O — XloTh- 2 uj 


OoO 


oo* 








* ClooD-Cl 2 


*0X0 Q_ ID U J • — i 


mZILlu 


o* 


o 


o 




* <>CDIT> O 


* — 2-j cl2_j 


U-OO ID 


CL* 




l~ 




* :^oo C) ID UJ 


* uj-*<t <r uj uj q_ 


LL ►— « ) 


* 


► — < 






* or D-JJ2 


* X* ZIi^JS 


LUh- X< 


-J* 




o 




* oujujujO or 


— * -Oh- — uj»— 


o<*-> 


<* 


< 


o 




* LLXtOQUD 


2* -DO 


ooo 


>* 








* h- o h- 


00 * 00 UJ<2^UJOO 


o— o 


<* — 








* < CL UJ LU UJ 


<* 


UJLU^UJ 


2* M 


2 


2 


♦> 


* or 5" x X or 


** oOllOO — 


X h- 


* 1 


oo 


oo 




* LOOOl — h- 


x* ok- a 


hlLO< 


•** o 


> 


> 




* 2X0 X 


•*-# CLoOoOOOlUXoOh- 


020 


2* - 


< 


o< 


< 


< or U'U LL CrT 


2* 5ThMh>-^ZD 


x c o 


~* X- 


• 


o # 




* G2TOC< 


-* OoO ^OCCi—o 


□a: coo 


H-4r l n^“ 


LUO 


OOUI 


< 


— * x CD 


<*o >- a:h-< I 


UJ _J 


2* <0 


-J< 


-J 


< 


2* or — 2 u iuoo 


— * UJ_J_J_Jh-<l-h- 


UJCDX <X 


<* w. ^ 


• ►— < 


a • 


1 


•• * x y— Uj 2 z or 


H-* Q<r2-JoO-£-2*~« 


2Z J(j 


O* OOO 


— o 


i — * 




2 * CL ID u <<X 


-J* tDO<UJ ODt 


<o< 


* * r-l 




— o 


-J 


r * h- cn 2 2 IC 


O* UJ*-H o . o 


2ZIU! 


00 * _J *-H 


2 O— 


Ooc 




CD* LU— 00 


— i* 2 q:di >- x 


X 


Uj ^ 


2 2 2 


O 


♦ 


** 2K-LOU 'a' 2 


* ^hZhoo<xo: 


lUUJlUl— 


_J# UJcoO 


► — * * — » *■ r 


* — • o ► 


-J 


M * — OO < X X < 


lu* h-xK^-.ar.— o 


xxx 


-J* CS < 


< r\j . — * — * c\i 


— — * O 


w 


* 1 — CO jChh - 


2*rx^:cD3: cCarL l 


y-y-y~<S) 


— * 2 * 


II - ii ii 


2<r^ n 


< 


>* ID ID X 


-<* Ox — <h- 


— . 


O* h-O 2 


►— i ►— • o <r iu o 


• — “0 ^ 


II UJ 


-j* aoo— cooox 


H- * a: >~ a: — <r <ro 


0 O 0 O 0 O 


* 1 ' ► — ' oO 


— oo ii ID 


r-H | OO - | 


— O 


oo* or < — — h- 


ID* cDooxi<^iu2lu 


►h > — * » — « 2 


X^UoO<hO<OcQ^ZO 


1 OQD— 


* COO — 


0*0 a- x o 


00 


CD* •”•2 II 1 


ono<“5 — oo . *r — o 4- 


rr — a:>-* iDor 


ar* oo o cl x<uj2 


<2Z< 


*— ILUZ2f<lOHQ ^h-CO + _jO<I(\JrH 


•* h- O CJC * oo < x < CD 


CO* LUID>-— 00X< 




O * CL 2 OO II 


< ^-2 


11 II J2hh^- 2 ►— 


ID* oOO <ctT 1 -cD 




UJ*M— .>5QmOu.^ oo II LL LL <f O — 1 


-pLL 


2* oo^ror 


oO* *— * 2 X! Cl 1— UJ 




O* hQ<2qDC^<UDJmm<qz 


< o ac UJ * M Uh 


* i<Xu:<iiLZ 




O* 








«• x<< 


-km-coi— eXf-co — 




O* 


O 




oo 


* K-CD2 


* 




* 


o 




oo 


* 


* 




* 


r-4 




CM CD 


* 


* 




* 











UOUUUUUOOUUUUUUUU oidooooo 



114 







OOOOOQOOOOOOOO^OOOOOO 
(MoOvtLnsor-coa>o— <rsj^Ov**irwor^-cocr' 0 -HrNj 
cvlfMosjfvjrvjrMfVjrvj^fOrnmr'^rn^fnrom^^^* 
O 0000 00 OO 000 00 oO 000 

CO to CO CO CO co CO co 00 CO co co CO CO CO CO uO CO CO CO co 

0 0000 no 0000 00a 000 0000 






-J 


_J -J _J -J -J -J _J -J -J -J O _J -J -J _J -J —1 _J -J -J 






* 






O 


->< 






LU 








or 








¥—i 


# 






3 


* 






O 


# 






LU 


# 






a: 


tt 








CM# 






co 


I s -* 






►— t 


O' # 








H 4 C 






a: 


* 






0 


>-# 






Y- 


<* 






0 


Xtt 






LU 


* 






> 


► * 








O* 




21 


O 


O* 




LU 


<r 


O* 




h- 


0 


x* 




co 


_j 






>- 2 T 




OO* 




COllJ 


or 


* 




h- 


UJ 


LLJ^ 




UJCO 


(X 


►- * 




X>- 




<* 




h-co 


> 


3 * 






3 


O* 




2 UJ 


CO 


<* 








<£* 




h- 


LL 


O# 




CO 


O 


H* 




Xu. 




co# 




00 


_J 


0 # 






-J 


CL.y. 




h-X 


<t 


* 




<n- 


O 




CL 


30 




<# ~* 


CL 


Q-< 


UJ 


>* O 


►-* 


LUS 


X 


< * 


CD 




0 


X# — ‘ 


* 


LLO 






O — 


OX 




CD 


O 0-5 


< 




X * - 


O O - 


or co 


_J 




Sf»— • O 


LU 


0 


# sf 


CO OCL 


COLL 


-J 


X* MO 


0 * 1 — ►— 


>T-J 




<* •* 


H- *— 


3 < 


0 


o> 0 


*0 —* 0^ 


XX 


k- 




O * - CO 1 






CO* r-J 


o*—* *— • 


LL LU 




LU* *-* 


N 

M 

( 

N 

( 

n 

) 

M 

) 

0 


XX 


-J 


3 * M 


► f — M * M * r — t r O O'- 


Y-Y- 


< 




— » CNJ X 1 ^^CNJ + CMCM •-> 




0 


X 


11 11 • - 11 — 11 a: 11 + tu*— 


cOcO 




O* O 


■— “5 Y— IIU» — « * — ' . — J «— f » — ' ~ 0 . icOlu 


M » — < 


LU 


* »-• 


— — 3 — +— 1 • II 3 




X 


CO 


^OmO 1 •a^XOcDO-JcDOcLO — XX 




0 


CO# x 


1 O-OOJ 11 *-*0 II 0 1 11 O^CLO— 'CL 


22 




*iu 


2: <}- co m- 11 — • cl h-o 




cL 




II II •— JZ — II cL 0 — ^ 2 H 0 




LU 


LU # — 


2 GhO II LL^OO*-OcL— 'OcLLL— OLL'X 




H 


O# QZOcOQ-J>~«cDUOcDO'-<CDOi--imcD(JCiI UJ 




LL 


O# 






< 


O* 


OO O 






* 


OO O 






* 


vt Ln no 


0000000 





115 



ocoooo°oooonoooooooooo r 'oo r )oooooo 

OO0000O0O >— < — < « — < • — c ■— < »-j < #— < * — < c\! cm cm cm c \j r\j cm c\j cj cm rn O r°i 
OOOO 00^0 0^00000 OOOOO 0 0^0 0^00 Oo^o 
DC cc dc cc ct.cc i€ a. oc o' a' u a cc a: cc cc or o' or cr' a. cZ a ft' o' ^ a: a: or on a: 
<<T<<<<<<<T<<< <<<<<<<<<< <<<X<<<<<T<< 

^ >1 :g tc :r x s; 2: >_ :• i'l'sxszssss 7 l 2: 2: r: 3 T ?- ?:srz2:z 





OO 






- •» 






71 *-H 






LUX 






LU *• 






OO 






DC - 












LUO" 






UJ ^ 






LL LU 






CL • 






LU CO 






LU O' 






ceze^ 






LL UJ 






ox 






LU CD 












CL>^ — 






0 x 






OX 






LU - 






IL 7 M 






5:0- 






O 






LUOUJ 






LU *- 






HZh 






^:o- 






LQ < 






UJOLU 






>-- 2 






1 - 21 - 






IS) + >-4 






l/> < 






XO 






>-► 7 






-JIAO' 






uo 






< -O 






xo 






0 X 0 






jina: 






•—XO 






< •CJ 






DC X 






oxu 






O **XI 






•^xo 






GT* - 






cc X 






•—re - 






UJ 






jOX 






X- - 






X-JCM 






CL - 






O-J— » 






c/jCJX 






O - 






.JC\J 






2 TU-- 






Z J • — c 






*— UJ 






—•O - 






LOl— 






IL . 






^Ul< 






2 UJ 






UJ|-^ 






U J uo l — 






><•— 






> LU < 






— 2: 0 






►-«) — z 






O*— •ct' 






o< — • 






00 






^O 






uoctro 






uo-^a: 


-»o 




UJOO 






iuOo 


II ^ 




xO 






K CLO 


_J « 


M 


<o>* 






<ou 


-M* 


— 1 O 


^ - 




XX 


7:0 


— CM 


rco 0 


~oc - 




00 


— o>- 


JC 


- - 00 


o<x 


CM CM 


;t;z 


0 ► 


• n 


^x a>o c\j 


CL —KM 


OO 


<< — 


OCDC ► 


O CL 


»• 1 *•*» x 0 0 ^ *~* 


OD-< 


0003 




o<x 




0 < rocz co OO 


Go - 


• • 


uo;zo 


OJ'M 


O * 


*C 0 O • 2 f IT 


OZG- 


• 4 r V 


OwZ 




— X 


^OOOl Ah- << 


- <UJ 


\\ 


OL 0 < 


- 0 • 


- lA 


cC -h- f-^ro N 


1 — 1 — 


» — « »— * 


— OO — 




O * 


< j^ino uozr 


x o< 


a Cl 


Mr r- uo 


X< LU 


LA 


^<-002. a. 0^ 


— tniuz: 


SI* * 


o<: r^o 


-ini-H 


7 " — * • 


ujum * -ft uooo ► ac »— 




-o< 


► O - 


iu ex'. OOO- - 0 


'-‘rocM 


<LOuOOO- IUZhOX 


2 T ^ LT\ CM n CNJ * # evo -0:0 1 ! - - 


— * ■* # A O CL ► -< 


II COO 


w hoh *•— cl or 


-- xoo 


*—<►—<* — 1 — * 


2: c<cl a: 


- 00 


h^^<MU ►— — t 11 II O^vOOvO 


1 


•— 11 n 11 0 — *0' 


0— 


DOn^ZO^ — 


— h- -0 


000 - 


00—* — —*. 


— h-xoo — H- ^ 


Om^ « • O O ►hH fNJ LU < X OO || |l ►— * 


Or-KMfOLU<vOOOUJ<Q^ 


DC — 1 LU 00 .'Z II — •> - 


hZSXh^NX- 


II - * - 


v - t : -oxi— 7:0 


CDCl X II — 0 O ft- - - 


O 0 0 2 ►— * 


^arx 


^OThD 


— u_o;z 11 -- 


cro- -oorz 11 


» 4 w w w 


c^oxxoa'OU’Z 


oowQaMa<o:uo 3 :ii.Luxoo<<Q:i r )ouo 


r< u_ •* ► 


O GC LL CL l U 




r^C\ f 






—<(M 




0 


O 


0 


0 


O 


OOO 


0 


O 


0 


0 


O 


OOO 


rH 


CM 


ro 


vf 


LA 


vOXCO 



116 



0000^)00^00 
r-< f\j ro ir\ vO r*- <*o o o 
ooooooo^o— < 
ooooooonoo 

OOOOOtAVloOCOOOoOLO 

cDcoajcocococCjcocnco 



00^00000000000000 

^rvjoOvMnor^ r oa'o»-™*f\iooNtu r >o^ 

00<^0000^>0rM— 4r— I 

OOOOOOO “*00^00^00^ 
I — > — I — t — I — h- I — b - I — I — H- h~ I hh Kh 
DDD^D*ODO)^^nO)DD:0 0) 

sx.rsyzs m ?; >. :r p >' :e :*_ p re 







o 






p 




f— ^ 






o 




f— < 






►— < 




w 






CD 




o 




O 


< 




•“ 




a: 




P 


«r-* 




► — < 


O 


w 


p o 






or-' 


H* 


^-1 




CD 


NN 


< 


P r^Jr-i 




o * 


CX- I 


O 


- 1 ~ 




O — — 


<ro 


_J 


uocn 




_j oct: 


•■ * 


— ll 






O 0~ 


ZI- 


— o 


COX — 




O O <N~ 


- i vt 




- 1 vf 




h- ■ CO 


<<>o 


••inro 


<<o 




- O-tt 


WW ► 


►Hr-lQ 


wW »■ 




O JC h-~ 


7COO 


— 1 o 


HCOO r-l 




O o 3 


UO* r-l 


<r ro # 


J "K r— < * 




a' n - 


P-Jr-I 


~Or-t 


3— !•— 1 ~ 




*-> *— « 0 1 — < 


cn< — 


7' 00 • g 


s<^z — 2 : 




O w — 


UJ< 


* CQ r-lO 


UJ< •> < ► 




• < — < 


tua: 


ex 


IUQC H * C\|-)UJ(\I+ P 4- 


2 PO II OCDN 


z p II II 


+ D + - •- 


— h-OO— + Z_'OL 


*-*!— O*— ' 1 1 


• — 5 r-t *— • UJ 


h- O 00 < < 




oO 


IS | ~ 1 Ow3 


DUO0 


• op ii ii p 


DUoOO | CD O 


*— < 


OPOO »OPP 


n^zoo<o^a: 


C^PO— II o II 


oi H ii 4* d: ii *—* ex 


a" -Jiu 


ii — i n a: 0 3 


q: _jujr\i ii — ex 1 




— —1 — — — 1— 3 


coa-s: 


CD Ol f^si *—h ) — O 


01 CL P |— •— 


o 


wO>- 11 '-MZhO 


32>- 


zozaaiuz 


Ds;«o^-oa:iuJ-ci:ii--aiu2 


oo— o<o<<<o:uj 


(/)^QO'“'UO 




» — < ' — 'O^^UOO^ LU 




o 






o o 




o 






o o 




r-t 






^ (M 



117 



SUBROUTINE STJF F ( A , B ♦ N , M) STF0010 



00 <~> 0 C >00 
Nr^sMAvO^-CO 
OO^OO^O 
00°000 o 

IL LL U_ LL L . U_ U- 
H H H h- H H 
Ul iS) UO t/l LO LO LO 



ooooo 

0 0r-<r\jf'0 

CD r— < r— I r— 4 r—4 

ooo^o 

U_ U_ IL U_ LL 



00 O O </) t/) v/) 



ooo^ooooooooo 
rHNrO'tin^NcOO'O^NfO 
OOOOOOOOO^r-i— «— i 

OOOOOOOo ooooo 

2 : i: 2 : z r' l: .-r ^ ^ 

CLCLCLO-aO-CLCLCLCLCLCLa. 



o 



< 

< 



7^ • * 



t— ICO 




X 




:s 




►— < 




r o 




a: 




<oc 




H 




z 




< 




IL 




Vi 




xo 












LL 




crct 




O 




H-ui 








<co 












o 




3 




V 




u.^ 




— 




o 








II 








OO 




v 




hZ 




3 




t: 




3 




uj • •* 




O 




HO 




o 




2LU 








o> 




o 




oo: 




IU 




UJ 




1— 1 




LUO0 




LL 




XUJ 




» — < 








o 




a. 


sO 


UJ 


— * 


UJ 




CL 


o 


o^ 


o 


OO 


r-k 


< — 


r — < «— «* 




H 


3 


Mr-1 O 


o 




CLL 


1 — iL 


H 


CO 


<X oO 


O co - 


2 : 


CO 


or 7^ 


• - ^ 


KH 


* 


-j v: 


X — 






hXD 


1 CO 


r 


m 


— 3 


<o co 


co 


► 


ora 


w r r 


CO 


o 


UJ 1 — o 


OOO < 


r 


r—t 


^<r 


*- —« <r 




•— 1 


— tzu. 


— J r— < * tar 


CL 


%-r » — 4 


H O 


< — 2 o > 


o 


<2 — 


3** 


UJ< * X * LL 


H 


< -CD 


or a: 


O' r-< * — LL 


O 


i — i CO 


OCcOUJ 


X H — 3 


IU 


71 II II 


cor co 


HO"- 11 < H 


> 


O — — 


3 vc 


— • 3 || O0 




— o 


00X3 


ooooo^z 


on 


uoOSZ 




mZOOo^> 


UJ 


Zo -CL 


O0nr 


juir^-i oq: 


o 


UJ(M — 3 


— H || 


CL2 : — HH 


< 


21 — hi 


x<r 


x:— oo— u 2 : 


3 


— 0<U 


1-22 


— ooocoa: iu 


CL 


QO<ari 




o 




o 




o 




o 




r~4 




C\J 



- I O 

ZO-fO 

► ► r— < O 

CDX-^ ~ 

- I CO — 

<< -o — >c 

" — ' - — ' " — ' > — < I 

UCOO • • 

D-I^O (\J 

Cl<w- 2^0 

UJ<^ * O - 
UJ cf C\j *— < O O 

2 2 • II CNi » — I O || iPH 

M hOh^ »mO m - — 
h — ■ — o r-'-'O r-w- 
oooo *o — h o*-h^ 

0hh2^OUJ<00lu<^ 

cc j id — — u— ix h nJ~ > — tod 

cQcl —a: —.arh- O 

D2^iL0^0U0a'0uuZ 
oo ►— O o 3 IL OQ IS IL CL UJ 

oo oooo 
oo oooo 

f~~< c\j ro^msO 



co 



— <\ J 

a 

► O . 

^oo 



ooooo OoO 



118 



LIST OF REFERENCES 



1. Arpaci , V. S., Conduction Heat Transfer , Addison-Weslcy 

Publishing Company, 1966. 

2. Berg, P. W. and McGregor, J. L. , E lementary Partial 

Differential Equations , Holden-Day, 1966. 

3. Crandall, S. H., Engineering Analysis, McGraw-Hill, 

1956. 

4. Haji-Sheikh, A. and Sparrow, E. M. , "Transient Heat 

Conduction in a Prolate Spheroidal Solid," Tran . 
ASHEHT , v. 88, p. 331-3, Augus L 1966. 

5. Holman, J. P., Heat Transfer , 2d ed . , Mc-Graw Hill, 1963. 

6. Irons, B. M. , "A Frontal Solution Program for Finite 

Element Analysis," I nternatio na l Journal for Numeri - 
cal Methods in Engineering , v. 2, p. 5-32, 1970. 

7. Johnson, D. E., "A Proof of the Stability of the Houbolt 

Method," AIAA Journal , v. 4, Mo. 7, p. 1450-1451, 
August 1966. 

8. Newmark, N. M. , "A Method of Computation for Structural 

Dynamics," ASCE , Journal of Engineering Meehan ics 
Division , v. 85, No. EM 3 , p. 67-64, July 1959. 

9. Stroud, A. H. and Secrest, D. , Gaussian Quadrature 

Formulas , Prentice-Hall Inc. , 1966. 

10. Zienkiewicz, 0. C., The Finite Element Method in 

Engineering Science , McGraw-Hill, 1971. 

11. Zienkiewicz, 0. C. and Parekh, C. J., "Transient Field 

Problems: Two-Dimensional and Three-Dimensional 

Analysis by Isoparametric Finite Elements," Inter - 
national Jourmil for Numerical Methods in Engineer - 
ing , v. 2, No. 1, p. 61-71, January-March 1970. 



12. Users Manual , W. R. Church Computer Center, U.S. Naval 

Postgraduate School, Monterey, California, 1970. 

13. Private communication from Gilles Cantin. 



119 



INITIAL DISTRIBUTION LIST 



No. 



1. Defense Documentation Center 
Cameron Station 
Alexandria, Virginia 22314 

2. Library, Code 0212 
Naval Postgraduate School 
Monterey, California 93940 

3. Professor Gilles Cantin, Code 59Ci 
Department of Mechanical Engineering 
Naval Postgraduate School 
Monterey, California 93940 

4. Asst Professor David Salinas, Code 59Zc 
Department of Mechanical Engineering 
Naval Pos tgx'aduate School 

Monterey, California 93940 

5. LCDR Girard T. Lev; 

Supervisor of Shi pbuilding. Conversion , 
and Repair, USN 

Pascagoula, Mississippi 39567 

6. Professor Paul F. Pucci, Code 59Pc 
Department of Mechanical Engineering 
Naval Postgraduate School 
Monterey, California 93940 

7. Professor Robert E. Nev/ton , Code 59Ne 
Department of Mechanical Engineering 
Naval Postgraduate School 
Monterey, California 93940 

8. Professor John E. Brock, Code 59Bc 
Department of Mechanical Engineering 
Naval Postgraduate School 
Monterey, California 93940 

9. Dr. Arnold Adicoff 
Naval Weapons Center 

China Lake, California 93555 

10. W. R. Church Computer Center, 0211 
Naval Postgraduate School 
Monterey, California 93940 



Copies 

2 

2 

5 

5 

3 

1 

1 

1 

5 

2 



120 



UNCLASSIFIED 





1 " ~ DOCUMENT CONTROL DATA - p. r 

1 /« t Sllic.lmn cl title, body ol eh, feci end indc,lr C eenclelion mull he • > 


t 

if trad when the ov»te ft report 1% ctrmltied) 


ting a C T l v i 1 r (t orporaie author) 

■val Postgraduate School 
'nterey, California 93940 


2 ,. REPORT SLCUR TV CLASSIFICATION 

Unclassified | 


2b. CROUP 


i title 

THREE DIMENSIONAL SOLUTION OF THE TRANSIENT FIELD PROBLEM 
ING ISOPARAMETRIC FINITE ELEMENTS 


p T)V f NOTES (T>r* of report end. InC tun ,ve datea) 

star's Thcs i s ; December 19 7 2 __ 


; n5 l (Fit H name, middle Initial, teat name) 

rard T. Lew; Lieutenant Commander, United States Navy 


\ DATE 

cember 1972 


7». TOTAL NO. OF PACES 7b. NO OF REFS 

3 22 13 


IACT OR GRANT NO 
CT NO 




««, ORIGINATOR** REPORT NUMBERS) 


6 b. OTHER REPORT NO(SI (Any other numb* t a that may be a a atoned 
thta report) 



II BU T I ON STATLmCNT 



bproved for public release; distribution unlimited. 



EMENURV NOTES 



12. SPONSORING MILITARY ACTIVITY 



Naval Postgraduate School 
Monterey, California 93940 



An isoparametric finite element formulation solving a 
Bneral form of the transient field equation is presented in 
tis thesis. The formulation is developed for fields in 
iree-dimensional Euclidean space. A FORTRAN IV computer 
cogram employing double precision arithmetic, compact storage 
achniques, and providing the option of several numerical inte- 
ration methods and types of finite elements is presentee . 
ata input may be in rectangular, cylindrical, or spherical 
oordinates and variations in time integration step size are 
ermitted. Time dependent boundary conditions are not con- 
idered . 

Comparisons of theoretical and computer solutions ior a 
ariety of test problems demonstrate close agreement. Instruc- 
ions for use of the program are included. 



> f rj473 

C)t -807-681 1 



(PAGE I ) 



121 



UNCLASSIFIED 



Security Cle«*ificstion 






\ 



i n 



\ aw "7 

- j ° 








\ 



Thes i s 141635 

lCC4 Lew 

C) 1 A three dimensional 

solution of the tran- 
sient field problem 
using isoparametric 
f i ni te elements. 



thesL604 

A three dimensional solution of the tran 




3 2768 002 11866 3 

DUDLEY KNOX LIBRARY 



